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Largazole  is  highly  selective  against  tumor  cells  but  not  normal  cells.  We  will 
develop  more  potent  and  selective  small  molecules  to  validate  the  concept  that  dual 
specificity  inhibitors  are  better  anti-cancer  drugs.  Our  studies  are  expected  to 
provide  novel  ideas  for  designing  more  effective  therapeutics  for  breast  cancer 
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1.  INTRODUCTION 

Histone  deacetylases  (HDACs)  and  histone  acetylases  (HATs)  are  key  players  in  regulating  transcription  and 
histone  homeostasis(8).  Transcription  of  tumor  suppressor  proteins  is  frequently  silenced  in  tumor  cells  due  the 
hyper-  or  aberrant  activity  of  HDACs.  Accordingly  inhibiting  histone  deacetylation  may  re-activate 
inappropriately  silenced  genes  and  may  be  able  to  “reverse”  malignant  changes(2).  Inhibitors  of  histone 
deacetylase  enzymes  (HDACi)  have  recently  attracted  substantial  attention  as  potential  anti-cancer  drugs.  The 
selective  degradation  of  many  regulatory  proteins  in  eukaryotic  cells  is  mediated  by  the  ubiquitin  system(5). 
Proteins  targeted  for  degradation  are  usually  covalently  ligated  to  a  polyubiquitin  chain  and  subsequently 
eliminated  by  the  26S  proteasome.  Ubiquitination  of  proteins  is  carried  out  by  a  multi-enzyme  complex 
consisting  of  El  (ubiquitin  activating  enzyme),  E2  (ubiquitin  conjugating  enzyme)  and  E3  (ubiquitin  ligase)  (5). 
The  final  product  of  this  reaction  is  fonnation  of  a  polyubiquitinated  protein  with  attachment  of  an 
ubiquitination  through  an  isopeptide  bond  to  an  epsilon-amino  group  of  certain  Lys  residues  in  the  interior  of 
the  substrate.  There  is  only  one  ubiquitin  El  enzyme,  more  than  fifty  ubiquitin  E2  and  perhaps  thousands  of 
E3  enzymes  in  human  genome.  E3  often  controls  the  specificity  and  timing  of  substrate  ubiquitination  (5). 
Both  HD  AC  inhibitors  and  ubiquitin-proteasome  inhibitors  have  found  applications  in  treating  specific  type  of 
human  tumors.  However,  either  type  of  inhibitor  alone  does  not  appear  to  exhibit  a  broad  spectrum  of 
inhibition  in  treating  a  variety  of  human  cancers.  These  observations  have  prompted  investigations  using  a 
combination  of  both  types  of  inhibitors  in  anti-tumor  studies.  It  was  found  that  bortezomib  killed  multiple 
myeloma  cells  more  efficiently  when  combined  with  histone  deacetylase  inhibitors(l  1).  Thus,  administering 
two  inhibitors  simultaneously  targeting  both  pathways  could  be  a  feasible  therapeutic  strategy  for  cancer 
treatment. 

The  goal  of  this  study  is  to  test  that  dual-specificity  small  molecules  capable  of  targeting  two  or  more  aberrant 
signaling  pathways  associated  with  human  cancers  will  be  more  efficacious  in  suppressing  human  tumors.  We 
found  that  Largazole,  a  cyclic  depsipeptide  natural  product  isolated  from  the  marine  cyanobacterium  by  Luesch 
and  coworkers  (15),  has  remarkable  potency  in  selectively  inhibiting  the  proliferation  of  breast  cancer  cells 
without  significant  effects  on  normal  breast  mammary  epithelial  cells.  We  found  Largazole  also  inhibits 
ubiquitin  El.  The  inhibitory  activity  of  these  small  molecules  on  ubiquitin  conjugation  has  been  traced  to  their 
inhibition  of  the  ubiquitin  El  enzyme.  To  further  dissect  the  mechanism  of  El  inhibition,  we  analyzed  the 
effects  of  these  inhibitors  on  each  of  the  two  steps  of  El  activation.  We  show  that  Largazole  and  its  derivatives 
specifically  inhibit  the  adenylation  step  of  the  El  reaction  while  having  no  effect  on  thioester  fonnation.  El 
inhibition  appears  to  be  specific  to  human  El  as  Largazole  ketone  fails  to  inhibit  the  activation  of  Ubalp,  a 
homolog  of  El  in  Schizosaccharomyces  pombe.  Moreover,  Largazole  analogs  do  not  significantly  inhibit 
SUMO  El.  Thus,  Largazole  and  select  analogues  are  novel  classes  of  ubiquitin  El  inhibitors  and  valuable  tools 
for  studying  ubiquitination  in  vitro.  This  class  of  compounds  could  be  further  developed  and  potentially  be  a 
useful  tool  in  cells.  Using  breast  cancer  3D  culture  model  and  xenograft  models,  we  show  that  Largazole  and 
our  newly  made  Paragazole  have  potent  antitumor  activity  toward  triple  negative  breast  cancer.  Our  ultimate 
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goal  is  to  demonstrate  that  dual  targeting  of  both  pathways  is  the  underlying  mechanism  for  the  potency  and 
selectivity  of  Largazole  for  breast  cancer  cells. 

2.  BODY — Studies  and  Results 

Three  specific  aims  were  proposed  in  the  original  application.  We  describe  our  progress  in  the  context  of 
approved  SOW. 

Aim  1.  Synthesize  derivatives  of  largazole  with  potentially  improved  molecular  properties  and  improved 
selectivity  for  transformed  vs  non-transformed  cells. 

Task  1.1  Synthesize  initial  round  of  largazole  analogs  for  SAR  testing  (Phillips,  Months  1-6) 
(Completed) 

Task  1.2  After  initial  SAR  data  is  obtained,  design  and  synthesize  a  second  round  of  focused  largazole 
analogs  (Phillips,  Months  7-12)  (Completed) 

Task  1 .3  Write  and  submit  manuscripts  describing  the  intial  phase  of  SAR  studies  as  well  as  annual  report 
to  CDMRP  (Liu  and  Phillips,  Month  12)  (Completed) 

Task  1.4  Synthesize  larger  quantities  of  selected  largazole  derivatives  for  detailed  testing  (Phillips, 
Months  14-16)  (ongoing) 

Task  1.5  Synthesize  largazole  derivatives  targeted  to  specific  cancer  cells  e.g.  folic  acid  derivatives 

(Phillips,  Months  17-23)  (ongoing) 

Task  1.6  Testing  the  activity  of  largazole  derivatives  using  HDAC1  enzymatic  assay,  p27  ubiquitination 
and  El  thiolester  assays  (Liu,  Completed). 

Results 

Largazole  stabilizes  GFP-p27  expression  in  Kipl6  cells.  A  hallmark  of  many  advanced  cancers  is  an 
excessive  degradation  of  the  cyclin-dependent  kinase  inhibitor  p27,  which  is  directed  by  SCFSkp2-mediated 
ubiquitination.  Hence,  stabilization  of  p27  degradation  represents  a  rational  approach  in  cancer  therapeutics. 
To  identify  small  molecule  inhibitors  that  can  stabilize  p27Kipl,  we  performed  a  screen  of  -3000 
compounds  from  NCI  DTP  diversity  set  along  with  several  natural  products  in  our  collection.  For  the  cell 
based  screen,  we  generated  a  mink  lung  epithelial  cell  line  (Kip  16)  stably  expressing  p27  that  was  cloned  in 
frame  with  green  fluorescent  protein  (GFP).  The  resulting  N-terminal  GFP-p27  fusion,  detectable  by 
fluorescence  microscopy,  was  used  to  determine  the  levels  of  p27  expression  upon  treatment  of  cells  with 
the  compound  libraries  in  96-well  fonnat.  Much  to  our  surprise,  the  most  potent  hit  that  emerged  from  this 
screen  was  the  natural  compound  Largazole  (Figure  1),  which  was  first  described  by  Luesch  and  coworkers 
(15)  and  subsequently  synthesized  in  several  laboratories  including  ours  (1,  3,  10,  14,  15,  18,  19).  Largazole 
induced  a  robust  and  highly  unifonn  upregulation  of  GFP-p27  at  concentrations  as  low  as  1  nM  (Figure  2a) 
as  compared  to  the  expression  levels  after  treatment  with  the  proteasome  inhibitor  MG  132.  We  did  not 
observe  an  increase  in  GFP-p27  expression  upon  treatment  with  the  vehicle  control  DMSO.  This  result 
suggests  that  Largazole  can  stabilize  GFP-p27  expression  in  cultured  cells. 

Largazole  and  select  analogues  inhibit  the  in  vitro  ubiquitination  of  p27  and  Trfl.  Initial  investigation 
into  the  mechanism  of  Largazole  indicated  that  the  compound  stabilized  the  expression  of  p27  in  cells. 
Since  the  concentration  of  cyclin-dependent  kinase  inhibitor  p27  is  mainly  regulated  at  the  protein  level  by 
increased  polyubiquitination  and  subsequent  proteasomal  degradation,  we  hypothesized  that  Largazole  and 
synthetic  analogues  stabilize  p27  by  inhibiting  the  ubiquitination  pathway  (9,  13).  One  of  the  downsides  of 
cell  based  assays  is  that  the  effects  observed  may  be  attributed  to  the  influence  of  multiple  pathways.  For 
example,  inhibiting  the  proteasome,  elevating  transcription  of  GFP-p27,  or  inhibiting  Cdk  activity  can  also 
lead  to  an  increase  in  p27  expression.  To  tease  out  the  mechanism  and  action  of  Largazole  on  p27 
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stabilization,  we  decided  to  test  the  effect  of  Largazole  on  p27  ubiquitination  in  a  fully  reconstituted  system 
in  vitro  (16,  17).  To  test  if  Largazole  affects  p27  ubiquitination  in  vitro,  we  incubated  Largazole  (L)  with 
p27,  ubiquitin  El,  E2,  SCFSkp2,  and  Cksl.  As  shown  in  Figure  1,  adding  Largazole  significantly  reduced 
polyubiquitinated  p27,  suggesting  that  Largazole  can  block  p27  ubiquitination.  Since  Largazole  is  known  to 
be  a  histone  deacetylase  inhibitor  and  has  a  thioester  moiety  that  links  an  aliphatic  chain  to  the  core,  we 
decided  to  test  whether  inhibition  of  p27  degradation  can  be  linked  to  its  histone  deacetylase  inhibitory 
activity.  The  structure-activity  relationship  for  Largazole  is  relatively  well  understood  (12).  Therefore  we 
next  tested  a  series  of  Largazole  analogues  to  get  a  preliminary  structure-activity  relationship  on  p27 
ubiquitination.  To  investigate  this,  Largazole  ester  (E),  Largazole  ketone  (K),  Largazole  macrocycle  (M), 
and  seco- Largazole  (S)  were  tested  in  an  in  vitro  p27  ubiquitination  assay  (Figure  2b).  We  also  added  the 
HD  AC  inhibitor  Trichostatin  A  (TSA),  the  structure  of  which  can  be  found  in  Figure  1,  to  the  assay  to 
determine  whether  or  not  other  HD  AC  inhibitors  affect  p27  ubiquitination.  We  observed  that  Largazole  (L), 
Largazole  ketone  (K),  and  Largazole  ester  (E)  inhibited  the  ligation  of  ubiquitin  onto  p27;  however,  the  M 
and  S  analogues  and  TSA  failed  to  inhibit  the  ubiquitination  of  p27  (Figure  1).  This  result  suggests  both  the 
macrocycle  and  aliphatic  chain  are  required  for  ubiquitin  El  inhibition.  Furthermore,  the  result  also  suggests 
that  the  thioester  moiety  of  Largazole  is  not  required  for  inhibition,  because  the  ketone  and  ester  analogues 
were  equally  potent  in  blocking  p27  ubiquitination.  In  addition,  El  inhibition  is  unrelated  to  HD  AC 
inhibitor  activity  of  Largazole  as  both  ketone  and  ester  failed  to  inhibit  HD  AC  and  TSA,  a  known  HD  AC 
inhibitor,  does  not  block  p27  ubiquitination  in  vitro.  Prior  to  ubiquitination,  p27  is  phosphorylated  by  the 
Cdk2-CyclinE  complex.  We  carried  out  an  in  vitro  p27  phosphorylation  assay  in  the  presence  of  either 
DMSO  or  Largazole  in  order  to  test  whether  or  not  the  decrease  in  p27  ubiquitination  was  due  to  the 
inhibition  of  the  Cdk2-CyclinE  complex.  We  observed  that  Largazole  does  not  inhibit  the  phosphorylation 
of  p27  (Figure  2b);  therefore,  the  inhibition  of  p27-ubiquitin  conjugation  is  due  to  an  inhibition  of  the 
ubiquitination  process  rather  than  phosphorylation  step. 
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Figure  1.  Chemical  structures  of  Largazole,  synthetic  analogues,  Figure  2.  Largazole  stabilizes  p27  expression  in  Kipl  6  cells  and  inhibits  p27 
and  Trichostatin  A  synthesized.  ubiquitination  in  vitro  but  not  phosphorylation  of  Cdk2-CyclinE. 


Even  though  there  is  a  limited  number  of  proteins  in  the  reconstituted  p27  ubiquitination  system  in  vitro, 
tracing  the  real  target  of  Largazole  is  still  quite  challenging.  Fortunately,  we  have  previously  established 
another  reconstituted  in  vitro  ubiquitination  assay  of  Trfl  with  SCFFbx4  (20).  There  are  a  few  overlapping 
components  between  these  two  assays.  The  effect  of  Largazole  on  Trfl  ubiquitination  should  offer  some 
insight  as  to  where  Largazole  might  target.  To  study  the  specificity  of  Largazole,  we  added  Largazole  ester 
to  an  in  vitro  Trfl  assay  and  found  that  Largazole  ester  inhibited  the  ligation  of  ubiquitin  onto  Trfl  in  a 
dose-dependent  fashion.  Since  Trfl  and  p27  require  different  E2  ubiquitin-conjugating  enzymes  and 
different  E3  ubiquitin-ligating  recognition  subunits  in  order  to  carry  out  each  ubiquitination,  we 
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hypothesized  that  Largazole  and  select  synthetic  analogues  inhibit  a  step  common  to  both  ubiquitination 
pathways. 


Largazole  ketone  inhibits  ubiquitin  El  activation.  In  vertebrates,  there  exists  only  one  known  ubiquitin- 
activating  El  enzyme,  UBA1.  Since  both  p27  and  Trfl  can  be  ubiquitinated  in  the  presence  of  UBA1,  we 
hypothesized  that  the  inhibitory  activity  of  Largazole  is  due  to  the  deactivation  of  El.  To  test  this 
hypothesis,  we  incubated  Largazole  and  Largazole  ketone  with  recombinant  El  prior  to  carrying  out  an  in 
vitro  thioester  assay  we  described  previously  (7).  The  presence  of  a  fluorescence  signal  in  the  thioester 
assay  suggests  the  fonnation  of  El -ubiquitin  adducts.  The  dose  dependent  decrease  in  fluorescence  indicates 
that  Largazole  and  Largazole  inhibit  the  fonnation  of  El -ubiquitin  adducts  (Fig  3  AC).  The  dose-response 
curves  generated  from  Figure  3BD  suggest  an  IC50  of  approximately  29  pM  and  25  pM,  respectively. 
Activated  ubiquitin  is  nonnally  transferred  to  ubiquitin  conjugating  enzymes  (E2).  If  El  activity  is  inhibited, 
we  expect  to  see  that  defects  in  El  activation  should  impair  the  attachment  of  ubiquitin  onto  Cdc34  (E2). 
To  further  validate  El  inhibition,  we  included  Cdc34,  the  E2  enzyme  required  for  p27  ubiquitinaton,  in  the 
El  reaction  mixture.  As  shown  in  Figure  3EF,  in  the  presence  of  ATP,  fluorescent  ubiquitin  is  transferred 
to  Cdc34  indicated  by  the  presence  of  a  fluorescent  Cdc34  band  on  the  gel.  Upon  incubation  with  E2, 
Largazole  or  Largazole  ester  reduce  the  amount  of  ubiquitin  molecules  that  are  transferred  from  El  to  E2  in 
a  dose-dependent  fashion.  This  result  is  consistent  with  the  notion  that  Largazole  or  Largazole  ester  inhibit 
El  activity. 
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Figure  3.  Largazole  (L)  and  largazole  ester  (E)  inhibit  ubiquitin  El  in  a  Figure  4.  Largazole  ketone  inhibits  the  adenylation  of  the  El 
dose  dependent  manner  in  vitro.  ubiquitin-activating  enzyme. 


A  potential  caveat  of  the  above  experiment  is  that  if  Largazole  or  Largazole  ester  also  blocks  the  transfer  of 
ubiquitin  from  El  to  E2  we  would  have  seen  the  same  result.  To  rule  out  this  possibility,  we  first  produced 
ubiquitin  charged  El  by  incubating  ATP  and  fluorescent  ubiquitin  for  30  min  at  room  temperature  followed 
by  the  addition  of  Cdc34,  which  was  also  contained  with  Largazole  or  Largazole  ester.  If  either  compound 
block  ubiquitin  transfer  from  El  to  E2,  we  would  observe  a  significant  decrease  in  Cdc34  fluorescence 
signal  regardless  of  the  order  of  compound  addition.  On  the  other  hand,  we  should  see  the  opposite  results. 
As  shown  in  Figure  3GH,  Cdc34  is  fully  conjugated  with  fluorescence  ubiquitin  when  Largazole  or 
Largazole  ester  was  added  after  generating  fluorescent  ubiquitin-El.  This  result  suggests  that  Largazole  or 
Largazole  ester  neither  blocks  the  transfer  of  activated  ubiquitin  from  El  to  E2  nor  promotes  hydrolysis  of 
ubiquitin  thioester. 

Largazole  ketone  inhibits  the  adenylation  step  of  El  activation.  El  forms  an  ubiquitin-adenylate 
intennediate  during  the  course  of  its  catalytic  cycle  (4).  Thus  the  mechanism  of  ubiquitin  El  activation  can 
be  studied  by  assaying  ATP:PPi  and  AMP:ATP  exchanges  (4).  Production  of  AMP  in  the  [a-32P]-AMP:[a- 

T9 

P]-ATP  exchange  assay  guarantees  that  a  thioester  bond  is  fonned  between  El  and  ubiquitin,  while  the 
release  of  PPi,  measured  by  the  [32P]-PPi:[y-32P]-ATP  exchange  assay,  signals  the  formation  of  ubiquitin 


adenylate.  To  further  dissect  the  mechanism  of  Largazole  inhibition,  two  nucleotide  exchange  assays  were 
carried  out  in  the  presence  of  Largazole  derivatives.  For  these  experiments  we  used  Largazole  ketone, 
which  is  similar  to  Largazole  and  Largazole  ester.  From  the  results  shown  in  Figure  4,  it  is  evident  that  the 
first  two  concentrations  of  Largazole  ketone  (100  and  50  pM)  inhibit  ubiquitination  of  El  similarly  and 
were  also  inhibitory  in  both  types  of  exchange  assays.  The  lack  of  a  [~“P]-PPi  signal  suggests  that  the 
adenylation  step  did  not  occur;  consequently,  ubiquitin  could  not  be  transferred  to  the  active  site  cysteine  to 
trigger  the  release  of  AMP.  Both  steps  of  the  El-catalyzed  reactions  can  be  measured  by  the  AMP:ATP 
exchange  assay.  The  lack  of  an  [a-32P]-AMP  signal  further  suggests  that  the  adenylation  step  is  inhibited  by 
Largazole  ketone.  Thus  Largazole  or  Largazole  derivatives  act  on  the  first  step  of  ubiquitin  activation 
pathway  by  blocking  the  fonnation  of  ubiquitin-adenylate. 


Selectivity  of  Largazole  ketone  against  SUMO  El  and  Ubalp. 

In  addition  to  ubiquitin,  there  exist  several  ubiquitin-like  proteins 
that  covalently  modify  other  proteins.  All  of  the  ubiquitin-like 
proteins  have  activation  pathways  similar  to  ubiquitin  (6).  In 
order  to  study  the  specificity  of  Largazole  to  the  ubiquitin 
pathway,  we  incubated  Largazole  ketone  with  SUMO-activating 
El  enzyme  prior  to  carrying  out  a  thioester  assay.  From  the 
results  in  Figure  5b,  we  found  that  Largazole  ketone  is  ineffective 
in  inhibiting  the  formation  of  El -SUMO  adducts.  From  the  dose- 
response  curve  generated  from  the  SUMO  El  fluorescence 
results,  the  IC5o  is  approximately  450  pM  as  opposed  to  10  pM 
for  ubiquitin  E.  Thus  Largazole  is  relatively  selective  in 
perturbing  ubiquitin  El  activation. 
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Figure  5.  Investigation  into  the  selectivity  of  Largazole 
ketone  (K).  A,  Largazole  ketone  (K)  fails  to  inhibit  the 
ligation  of  ubiquitin  onto  Ubalp,  a  homologue  of  UBA1 
from  S.  pombe.  Formation  of  Ubalp-ubiquitin  adducts 
was  determined  by  thioester  assay  utilizing 
fluorescein-ubiquitin. 


Ubiquitin  and  the  ubiquitin  El  enzyme  are  highly  conserved 
among  eukaryotes  (6).  Sequence  analysis  shows  a  45%  homology  between  the  human  ubiquitin-activating 
enzyme  El  (UBA1)  and  S.  pombe  El  (ptr3/Ubalp)  at  the  amino  acid  sequence  level.  To  test  whether 
Largazole  ketone  inhibits  the  S.  pombe  El,  we  carried  out  a  thioester  assay  using  Largazole  ketone  and  the 
ubiquitin  El  homologue  in  S.  pombe,  Ubalp.  The  results  in  Figure  5a  suggest  that  Largazole  ketone  fails  to 
inhibit  the  formation  of  El -ubiquitin  adducts  at  concentrations  less  than  1  mM.  Taken  together,  these  results 
suggest  that  Largazole  and  its  derivative  are  highly  selective  in  inhibiting  the  ubiquitin  El  enzyme. 


The  results  obtained  in  Aim  1  have  been  published  in  PLoSONE  Journal  earlier  this  year. 


Aim  2.  To  determine  whether  dual  inhibition  of  both  HDAC  and  ubiquitin  conjugation  is  responsible  for  the 
electivity  of  largazole  against  breast  cancer  cells  and  determine  which  HDAC  isoforms  render  breast 
epithelial  cells  sensitive  to  largazole. 


We  have  started  the  effort  of  cloning  of  all  HDAC  enzymes.  So  far  we  have  cloned  HDAC1,  2,  3,  4,  6,  8 
and  11.  We  will  continue  to  procure  clones  from  ATCC  to  complete  the  HDAC  enzyme  set  cloning  project. 
The  TGI  assay  on  MDA-MB23 1  and  a  battery  of  breast  cancer  cell  lines  have  been  ongoing. 


Task  2. 1  Cloning  of  17  HDAC  enzymes  into  lentivirual  vector  (Liu,  Months  1-4)  (Completed). 

Task  2.2  Construct  human  mammary  epithelial  cell  lines  (HME)  expressing  each  individual  HDAC 
enzyme  (Liu,  Months  5-16)  (Partially  completed). 

Task  2.3  Perform  cell  based  growth  inhibition  assays  using  MDA-MB231  and  HME  cells  using  largazole 
analogs  generated  in  Aim  1  (Liu,  Months  6-18)  (ongoing). 

Task  2.4  Measure  cell  penneability  of  lagarzole  analogs  using  parallel  artificial  membrane  permeability 

assay  (Liu,  Months  6-18)  (completed). 
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Task  2.5  Measuring  the  largazole  sensitivity  of  HME  cell  lines  expressing  HD  AC  enzymes  (Liu,  Months 
12-18)  (Completed). 

Task  2.6  Perfonn  siRNA  and  shRNA  experiments  for  informative  HD  AC  enzymes  in  HME  and  MDA- 
MB231  cells  (Liu,  Months  19-24)  (completed). 

Task  2.7  Write  the  manuscript  describing  dual  targeting  activity  of  largazole  (Liu  and  Phillips,  Month 
12,  completed). 


We  tested  a  panel  of  18  breast  cancer  cell  lines  from  a  heterogeneous 
group  of  breast  cancer  cell  lines.  In  our  initial  exploratory  studies,  we 
made  an  interesting  observation  that  triple  negative  breast  cancer 
(TNBC)  cell  lines  are  particularly  sensitive  to  Paragazole  inhibition. 
Since  TNBC  are  associated  with  a  shorter  median  time  to  relapse  and 
death  and  significant  unmet  medical  need  due  to  the  fact  that  these 
cancers  do  not  respond  to  endocrine  therapy  or  other  available  targeted 
agents,  we  decide  to  focus  on  testing  TNBC  cell  lines  instead  of 
broader  spectrum  of  breast  cancer  cell  lines.  We  used  a  panel  of  19 
breast  cancer  cell  lines  to  assess  the  proliferative  response  to  increasing 
concentrations  of  Paragazole  using  an  SRB  assay.  As  shown  in  Figure 
6,  a  majority  of  these  cell  lines  are  quite  sensitive  (<IC50=500  nM)  to 
Paragazole  treatment  although  some  are  more  sensitive  than  others. 
Several  TNBC  lines  are  inhibited  at  subnanomolar  concentrations  of 
Paragazole.  These  preliminary  data  suggests  that  Paragazole  might  be  a 
novel  agent  to  combat  TNBC. 


TNBC  vs.  Paragazole 
72h  SRB  proliferation  assay 


Figure  6.Growth  Inhibition  assays  of 
Paragazole  in  TNBC  lines. 


Next  we  picked  a  sensitive  cell  line  CAL- 120 
and  a  relatively  growth  resistance  cell  line 
Hs578t  to  stably  over  express  HD  AC  1,3, 6, 8  in 
them  using  lentiviral  mediate  gene  transfer.  Of 
four  enzymes,  only  overexpression  of  HD  AC  1 
has  observable  differences  in  Largazole  or 
Paragazole  response.  These  differences  are 
quite  small  and  more  rigorous  studies  are  needed 
to  establish  that  the  effects  are  reproducible  and 
statistically  significant  (Figure  7).  Preliminary 
data  does  appear  to  support  the  hypothesis  that 
overexpression  of  certain  HDAC  isofonn  may 
render  cells  more  sensitivity  to  Largazole 
inhibition. 
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Figure  7.  Growth  Inhibition  assays  of  Largazole  sensitivity  in  cell  lines 
overexpression  of  HDAC  1. 


Aim  3.  To  determine  the  chemotherapeutic  efficacy  of  largazole  to  inhibit  breast  cancer  growth  and 
metastasis  in  mice 

Task  3.1  Large  scale  synthesis  for  largazole  for  initial  testing  along  with  control  (Taxol)  (Phillips, 
Months  Completed) 

Task  3.2  Large  scale  synthesis  of  the  improved  lead  for  animal  testing  (Phillips,  Completed) 

Task  3.3  Determining  the  chemotherapeutic  effectiveness  of  largazole  to  prevent  mammary  tumor  growth, 
invasion,  and  angiogenesis  using  the  xenograft  animal  model  (Schiemann,  Completed) 

Task  3.4  Assessing  the  effects  of  largazole  to  prevent  mammary  tumor  metastasis  (Schiemann, 
Completed) 
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Task  3.5 


Testing  newly  improved  largazole  analogs  in  growth,  invasion,  agiogenesis  and  tumor  metastasis 
in  the  xenograft  animal  model  (Schiemann,  Months  12-23) 

Task  3.6  Write  and  submit  manuscripts  for  publication  and  final  report  to  CDMRP  (Liu,  Phillips  and 
Schiemann,  Month  24) 


Because  Dr.  William  Schiemann,  the  collaborator  on  this  project  has  moved  his  laboratory  to  Case  Western 
Reserve  University  in  Cleveland,  it  took  a  while  for  the  BRCP  to  complete  issuing  the  funding  for  this  aim 
to  us.  Because  of  his  move  we  have  to  rework  the  subcontract  and  obtain  approval  for  IACUC  and  animal 
protocols.  In  addition,  Professor  Schiemann  had  to  rebuild  his  research  team  at  Case  Western  Reserve 
University  since  many  experienced  researchers  in  his  lab  did  not  follow  him  to  Cleveland.  Despite  of  these 
setbacks,  Professor  Schiemann  group  managed  to  make  progress  on  this  project.  We  made  a  derivative  of 
Largazole  called  Paragazole.  In  the  last  report,  we  mentioned  that  Professor  Schiemann  had  obtained 
preliminary  results  showing  that  1)  Paragazole  Potently  Inhibits  the  Growth  of  Human  and  Murine  TNBC 
Organoids  (Figure  8);  2)  Paragazole  Potently  Inhibits  the  Growth  of  Human  and  Murine  TNBC  Organoids 
(Figure  8).  We  requested  a  no  cost  extension  to  allow  us  to  complete  the  remaining  SOW.  This  includes 
repeat  in  vivo  tumor  suppression  experiments.  During  the  no  cost  extension  period,  we  have  shown  that 
Paragazole  (10  mg/kg)  was  highly  effective  at  inhibiting  the  development  and  metastasis  of  human  MDA- 
MB-231  tumors  in  nude  mice.  Importantly,  the  lymph  node  metastasis  ofMDA-MB-231  cells  was  completely 
absent  in  Paragazole  treated  animals  (Figure  9). 


Fig.  8.  Paragazole  potently  inhibits  the  growth  of  human  and 
murine  TNBC  organoids.  Human  MDA-MB-231  (A)  and 
murine  4T1  (B)  TNBC  cells  were  incubated  in  the  absence 
(diluent)  or  presence  of  increasing  concentrations  of  Paragazole 
(0-2  mM)  over  a  span  of  1 1  days.  Differences  in  organoid 
growth  were  measured  longitudinally  by  bioluminescence  (top), 
or  by  bright-field  microscopy  (bottom).  Data  are  the  mean  (+/-) 
STD  of  organoids  treated  with  250  nM  of  Paragazole. 


Fig.  9.  Paragazole  inhibits  human  TNBC  tumor  growth  and 
metastasis  in  nude  mice.  Human  MDA-MB-231  cells  (2 
million  cells/injection)  were  engrafted  onto  the  mammary  fat 
pads  of  nu/nu  mice  and  allowed  to  develop  until  the  tumors 
reached  a  size  of  400  mm3,  at  which  point  cohorts  were 
treated  with  diluent  (NS),  Paragazole  (10  mg/kg, 
Paragazole).  (A)  tumor  growth  curve  (B)  size  of  tumors  (C) 
Weight  of  tumor  (D)  imaging  of  growing  tumor  in  animals 
(E)  metastasis. 


Overall,  our  studies  clearly  demonstrated  these  Largazole  and  Paragazole  are  very  promising  in  treated 
triple  negative  breast  cancers  in  vitro  and  in  vivo.  The  manuscript  describing  this  work  is  current  in 
preparation. 

The  Development  of  a  novel  high  throughput  computational  tool  for  studying  individual  and  collective 
cellular  migration 

Understanding  how  cells  migrate  individually  and  collectively  during  development  and  cancer  metastasis 
can  be  significantly  aided  by  a  computation  tool  to  accurately  measure  not  only  cellular  migration  speed, 
but  also  migration  direction  and  changes  in  migration  direction  in  a  temporal  and  spatial  manner.  We  have 
developed  such  a  tool  for  cell  migration  researchers,  named  Pathfinder,  which  is  capable  of  simultaneously 
measuring  the  migration  speed,  migration  direction,  and  changes  in  migration  directions  of  thousands  of 
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cells  both  instan  taneously  and  over  long  periods  of  time  from  fluorescence  microscopy  data.  Additionally, 
we  demonstrate  how  the  Pathfinder  software  can  be  used  to  quantify  collective  cell  migration.  The  novel 
capability  of  the  Pathfinder  software  to  measure  the  changes  in  migration  direction  of  large  populations  of 
cells  in  a  spatiotemporal  manner  will  aid  cellular  migration  research  by  providing  a  robust  method  for 
determining  the  mechanisms  of  cellular  guidance  during  individual  and  collective  cell  migration. 

3.  KEY  RESEARCH  ACCOMPLISHMENTS 

•  We  have  completed  synthesis  initial  round  of  largazole  analogs  for  SAR  testing. 

•  We  showed  that  Largazole  and  its  analogs  selectively  inhibit  ubiquitin  El  enzyme  activity  in 
vitro 

•  we  demonstrated  that  inhibitory  activity  of  Largazole  is  independent  of  its  inhibitory  activity 
towards  the  histone  deacetylase  enzymes 

•  Structure-activity  relationship  analysis  shows  that  the  thioester  bond  is  not  required  for  inhibition 
but  the  macrocycle  core  and  aliphatic  tail  are  required. 

•  Largazole  blocks  ubiquitin  activation  at  the  adenylation  step  and  without  perturbing  ubiquitin 
transfer  from  El  to  E2. 

•  We  show  that  Largazole  inhibition  of  El  is  highly  selective  as  it  does  not  inhibit  a  highly  related 
ubiquitin  El  enzyme  from  S.  pombe  and  is  almost  twenty  fold  less  effective  in  inhibiting  the 
activation  of  SUMO  E 1 . 

•  We  show  Largazole  represents  a  new  class  of  ubiquitin  El  inhibitor. 

•  We  show  that  Largazole  and  Paragazole  potently  inhibit  the  growth  of  human  and  murine  triple 
negative  breast  cancer  organoids. 

•  We  show  Largazole  and  Paragazole  exhibit  cytotoxicity  against  Human  MDA-MB-231  tumors  in 
mice. 

•  We  show  that  Pargazole  inhibits  triple  negative  breast  cancer  cell  line  Mda-mb231  metastasis  to 
lymph  nodes. 

•  Overexpression  of  single  HDAC  enzyme  isoform  is  insufficient  to  convert  sensitive  cells  to 
resistance  or  vice  versa. 

•  Development  of  a  novel  high  throughput  computational  tool  (PathFinder)  for  studying  individual 
and  collective  cellular  migration. 

•  Using  PathFinder  we  show  that  Pargazole  can  inhibit  MDA-MB21  cell  migration  in  response  to 
TGF-beta. 


4.  REPORTABLE  OUTCOMES 


Ungermannova  D,  Parker  SJ,  Nasveschuk  CG,  Wang  W,  Quade  B,  Zhang  G,  Kuchta  RD,  Phillips  AJ,  Liu 
X.  Largazole  and  its  derivatives  selectively  inhibit  ubiquitin  activating  enzyme  (el).  PLoS  One. 
2012;7(l):e29208. 

Chapnick  DA,  Jacobsen  J,  and  Liu,  X.  The  Development  of  a  novel  high  throughput  computational  tool  for 
studying  individual  and  collective  cellular  migration.  PLoS  One.  Submitted,  reviewed  and  acceptable  for 
publication  pending  minor  revision. 
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5.  CONCLUSIONS 


We  have  made  significant  progress  in  our  proposed  studies.  We  demonstrated  that  Largazole  is  a  new  class 
of  ubiquitin  El  inhibitor  and  the  activity  of  El  inhibition  is  independent  of  its  inhibitory  activity  toward 
HDAC.  It  is  possible  to  design  novel  dual  inhibitors  toward  both  pathways.  We  have  published  one  and  an 
additional  paper  has  been  submitted  and  requires  minor  revision.  We  expect  to  publish  at  least  two 
additional  manuscripts  on  this  project.  We  will  continue  to  pursue  the  goals  outlined  in  the  original 
proposal  in  hope  to  develop  a  more  effective  anti-breast  cancer  drug  candidate. 
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Abstract 

Protein  ubiquitination  plays  an  important  role  in  the  regulation  of  almost  every  aspect  of  eukaryotic  cellular  function; 
therefore,  its  destabilization  is  often  observed  in  most  human  diseases  and  cancers.  Consequently,  developing  inhibitors  of 
the  ubiquitination  system  for  the  treatment  of  cancer  has  been  a  recent  area  of  interest.  Currently,  only  a  few  classes  of 
compounds  have  been  discovered  to  inhibit  the  ubiquitin-activating  enzyme  (El)  and  only  one  class  is  relatively  selective  in 
El  inhibition  in  cells.  We  now  report  that  Largazole  and  its  ester  and  ketone  analogs  selectively  inhibit  ubiquitin  conjugation 
to  p27Klp1  and  TRF1  in  vitro.  The  inhibitory  activity  of  these  small  molecules  on  ubiquitin  conjugation  has  been  traced  to 
their  inhibition  of  the  ubiquitin  El  enzyme.  To  further  dissect  the  mechanism  of  El  inhibition,  we  analyzed  the  effects  of 
these  inhibitors  on  each  of  the  two  steps  of  El  activation.  We  show  that  Largazole  and  its  derivatives  specifically  inhibit  the 
adenylation  step  of  the  El  reaction  while  having  no  effect  on  thioester  bond  formation  between  ubiquitin  and  El.  El 
inhibition  appears  to  be  specific  to  human  El  as  Largazole  ketone  fails  to  inhibit  the  activation  of  Ubal  p,  a  homolog  of  El  in 
Schizosaccharomyces  pombe.  Moreover,  Largazole  analogs  do  not  significantly  inhibit  SUMO  El.  Thus,  Largazole  and  select 
analogs  are  a  novel  class  of  ubiquitin  El  inhibitors  and  valuable  tools  for  studying  ubiquitination  in  vitro.  This  class  of 
compounds  could  be  further  developed  and  potentially  be  a  useful  tool  in  cells. 
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Introduction 

In  humans,  protein  ubiquitination  is  a  dynamic  process, 
depending  on  a  tightly  regulated  balance  between  the  activity  of 
two  ubiquitin-activating  enzymes  (Els),  approximately  40  ubiqui- 
tin-conjugating  enzymes  (E2s),  and  hundreds  of  ubiquitin  ligases 
(E3s).  Protein  ubiquitination  and  subsequent  degradation  regulates 
almost  every  aspect  of  eukaryotic  cellular  function  including  cell 
cycle  regulation,  endocytosis,  signal  transduction,  apoptosis,  DNA 
damage  repair,  transcriptional  regulation,  and  many  others  [1]. 
Hershko  and  coworkers  discovered  that  ubiquitin  covalently 
modified  proteins  prior  to  their  degradation  in  rabbit  reticulocyte 
lysates  and  characterized  the  reaction  mechanism  [2].  They  first 
described  the  ubiquitin-activating  enzyme,  El,  that  carries  out  the 
ATP-dependent  activation  of  the  C-terminal  glycine  residue  of 
ubiquitin  prior  to  ligation.  In  the  first  step  of  the  El  activation,  the 
enzyme  forms  a  complex  with  ubiquitin  and  ATP  and  catalyzes 
the  adenylation  of  ubiquitin  and  successive  release  of  pyrophos¬ 
phate  (PPi).  During  the  second  step,  a  thioester  bond  is  formed 
between  the  C-terminus  of  ubiquitin  and  El,  subsequently 
releasing  adenosine  monophosphate  (AMP).  In  the  final  step  of 
El  activation,  additional  ATP  and  ubiquitin  are  recruited  to  the 
adenylation  site,  generating  a  fully  loaded  El  carrying  two 


molecules  of  ubiquitin.  The  activated  ubiquitin  is  then  transferred 
to  a  cysteine  in  the  active  site  of  ubiquitin  carrier  protein  E2,  also 
via  thiol  ester  linkage.  Some  E2  enzymes  transfer  ubiquitin  to 
acceptor  proteins  directly,  whereas  other  E2s  require  additional 
substrate  binding  proteins  known  as  ubiquitin  ligases  or  E3s 
[3,4]  .Through  this  mechanism,  ubiquitin  is  attached  to  proteins  by 
isopeptide  linkages  between  the  C-terminal  Gly76  of  ubiquitin  and 
the  s-amino  groups  of  lysine  residues  present  in  substrate  proteins. 
In  addition,  linkages  between  Lys48  of  one  ubiquitin  and  the  C- 
terminal  Gly  7 6  of  another  ubiquitin  ultimately  form  polyubiquitin 
chains  [5].  Once  polyubiquitinated,  proteins  are  targeted  by  the 
26S  proteasome  for  degradation. 

In  many  human  cancers,  the  ubiquitination  system  is  often 
destabilized.  For  example,  the  cyclin-dependent  kinase  inhibitor 
p27  is  mainly  regulated  at  the  protein  level  and  is  excessively 
degraded  in  approximately  50%  of  all  human  cancers  [6,7], 
Furthermore,  expression  of  p27  is  primarily  controlled  by 
polyubiquitination  via  the  SCFskp2  E3  ubiquitin  ligase  and 
subsequent  proteasomal  degradation  [8],  The  SCFskp“  is  a 
cullin-RING  ligase  (CRL),  which  is  comprised  of  RING-box 
protein  I  (Rbxl),  scaffold  protein  Cull,  linker  protein  Skpl,  and  F- 
box  protein  Skp2  [9].  In  order  for  the  ligase  to  function,  Cull  must 
first  be  covalently  modified  by  NEDD8,  an  ubiquitin-like  protein 
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[10-12].  Therefore,  an  observed  stabilization  ofp27  in  cells  could 
result  from  decreased  polyubiquitination  by  inhibiting  the 
neddylation  of  Cull  or  one  of  the  enzymes  required  for 
ubiquitination. 

Given  that  ubiquitination  influences  many  cellular  functions, 
malfunctions  in  the  pathway  play  a  role  in  the  pathogenesis  of 
human  neurodegenerative  disorders  such  as  Parkinson’s,  Alzhei¬ 
mer’s  and  Huntington’s  diseases,  as  well  as  cancer.  Inhibiting 
components  of  the  ubiquitination  system  seems  to  be  an  avenue  of 
therapeutic  development  with  clinical  applications  [13,14],  For 
example,  each  E3  ligase  targets  a  small  number  of  proteins  for 
ubiquitination,  which  makes  it  a  potential  target  for  highly  specific 
inhibitors  that  have  few  side  effects.  There  has,  however,  been  little 
success  in  developing  inhibitors  of  specific  E3  ligases  until  recently 
[14,15].  Also,  proteasome-inhibiting  compounds  have  been  a 
target  of  interest  and  were  originally  developed  as  tools  for  probing 
its  proteolytic  function  [16,17];  however,  these  inhibitors  were 
considered  as  possible  cancer  therapeutics  after  it  was  observed 
that  diey  induced  apoptosis  in  leukemic  cell  lines  [18-20], 
Although  inhibiting  the  proteasome  would  nonspecifically  inhibit 
the  entire  ubiquitination  system,  the  proteasome  inhibitor 
Bortezomib  has  fared  surprisingly  well  in  clinical  trials  and  is 
now  FDA  approved  for  the  treatment  of  relapsed  and  refractory 
myeloma  and  mantle  cell  lymphoma  [19].  Therefore,  inhibitory 
compounds  of  the  ubiquitin  system,  whether  they  are  specific  or 
nonspecific,  have  the  potential  to  be  important  therapeutics  for  the 
treatment  of  cancer. 

In  January  2008,  the  Leuscli  group  at  the  University  of  Florida 
identified  a  natural  product  they  named  Largazole,  which  was 
isolated  front  cyanobacteria  of  the  Symploca  genus.  They  examined 
the  compound  for  cytoxicity  against  cancer  cells  and  observed 
remarkable  antiproliferative  activity  in  transformed  mammary 
epithelial  cells.  In  addition,  they  showed  that  Largazole  preferen¬ 
tially  targets  cancer  cells  over  normal  cells,  which  makes  this 
marine  substance  an  important  synthetic  target  as  well  as  a 
potentially  valuable  cancer  chemotherapeutic.  Remarkably,  the 
structure  consists  of  several  unusual  features,  such  as  a  16- 
membered  macrocycle  containing  a  4-mediylthiazoline  fused  to  a 
thiazole  ring  and  an  octanoic  thioester  side  chain,  a  unit  rarely 
found  in  natural  products.  [2 1] .  Also,  Leusch  and  co-workers  first 
reported  the  total  synthesis  of  Largazole  and  determined  that  the 
molecular  basis  for  its  anticancer  activity  is  HDAC  inhibition 
[21,24],  Numerous  analogs  of  Largazole  have  been  generated  in 
efforts  to  understand  the  structure-activity  relationship,  and  it  has 
been  determined  that  the  thioester  moiety  is  required  for  HDAC 
inhibition  [21-32].  Here,  we  report  in  vitro  mechanistic  studies  that 
reveal  a  potential  role  of  Largazole  as  an  antagonist  of  the 
ubiquitin-activating  enzyme  El.  In  contrast  to  HDAC  inhibition, 
ketone  and  ester  analogs  of  Largazole  can  actively  block  the 
ligation  of  ubiquitin  onto  El,  indicating  a  differential  mode  of 
inhibitory  activity  since  the  formation  of  a  thiol  metabolite  is 
indispensible  for  El  inhibition.  More  explicitly,  Largazole’s 
presence  negatively  affected  the  formation  of  ubiquitin  adenylate, 
which  we  monitored  through  nucleotide  exchange  assay. 

Materials  and  Methods 

Construction  of  Kipl  6,  a  GFP-p27  Expressing  Cell  Line 

Mink  lung  epithelial  cells  expressing  GFP-p27  were  generated 
by  retroviral-mediated  gene  transfer.  pBabe-GFP-p27  amphotro- 
pic  virus  was  made  by  cotransfecting  pBabe-GFP-p27-Puro  with 
pCL-Ampho  in  293T  cells.  Viral  supernatant  was  collected  and 
used  to  infect  mink  lung  epithelial  cell  line  Mv  1  Lu  (CCL-64)  from 
ATCC  in  the  presence  of  8  pg/ml  polybrene.  Puromycin  was 


added  at  5  |tg/ ml  and  stable  clones  were  selected.  Each  clone  was 
subcultured  and  tested  for  GFP-p27  expression  in  die  presence  or 
absence  of  10  pM  MCI 32  (Calbiochem,  Darmstadt,  Germany) 
for  24  hours.  Clones  expressing  high  levels  of  GFP  in  the  presence 
of  MG  132  but  low  or  unndetectable  GFP  in  its  absence  were 
expanded.  Immunoblotting  using  an  anti-p27  antibody  (Santa 
Cruz  Biotechnology,  Santa  Cruz,  CA)  was  used  to  confirm  the 
expression  of  the  GFP-p27  fusion  protein  and  stabilization  of  GFP- 
p27  upon  MG132  treatment.  One  of  the  clones  used  for  all 
subsequent  studies  was  named  Kip  16. 

Largazole  Treatment  of  Kip  16  cells 

T otal  synthesis  of  Largazole  and  Largazole  analogs  is  described 
in  [24]  within  the  supporting  information  (including  copies  of 
spectra  of  all  compounds)  and  is  available  at  http://pubs.acs.org. 
Kip  16  cells  were  seeded  into  96-well  flat  clear-bottomed  plates  at 
40,000  cells/well  in  100  pi  medium  and  incubated  overnight  at 
37°C  in  a  humidified  5%  C02  atmosphere.  Largazole  was  then 
added  to  final  concentrations  ranging  from  1  pM  to  1  nM  in 
300  pi  of  fresh  medium.  0.3%  DMSO  and  1  pM  ofMG132  were 
used  as  negative  and  positive  controls,  respectively.  After  24  hours 
of  incubation,  the  medium  was  removed,  the  cells  were  washed 
twice  with  phosphate-buffered  saline  (PBS),  and  the  cells  were 
fixed  with  4%  paraformaldehyde  in  PBS  for  15  minutes  and  stored 
at  4°C  for  microscopy  evaluation.  Cells  were  visualized  with  a 
GFP  filter  set  using  a  10  x  objective  on  an  Eclipse  TE2000-S 
(Nikon,  Melville,  NY)  equipped  with  a  Photometries  camera 
(Roper  Scientific,  Tucson,  AZ). 

UBA1  and  His-cdc34  Purification 

Human  ubiquitin  El  (UBA1)  was  expressed  with  an  N-terminal 
GST  tag  fusion  by  means  of  recombinant  baculovirus  expression 
in  Hi5  insect  cells  using  the  pFastBacHTA  vector  (Invitrogen, 
Carlsbad,  CA).  The  cells  were  lysed  by  sonication  in  the  presence 
of  protease  inhibitors  in  a  buffer  containing  200  mM  NaCl, 
50  mM  Tris-HCl  pH  7.5,  1%  NP40,  1  mM  DTT,  and  1  mM 
EDTA.  Cleared  lysate  was  incubated  with  glutathione  beads 
(Amersham,  Sweden)  for  one  hour  at  4°C.  After  three  washes  with 
lysis  buffer,  untagged  El  was  produced  by  thrombin  cleavage.  The 
protein  solution  was  passed  through  a  S200  gel  filtration  column 
(Amersham,  Sweden),  and  UBA1  concentration  and  purity  was 
evaluated  by  SDS-PAGE  and  Coomassie  Blue  gel  staining.  The 
purity  was  generally  greater  than  90%  and  purified  UBA1  was 
aliquoted  and  stored  at  —  80°C  after  quick  freezing  in  liquid 
nitrogen.  N-terminal  hexahistidine  (His)-tagged  human  Cdc34  was 
cloned  into  the  pQE-30  vector  (Qiagen,  Valencia,  CA)  and 
expressed  in  Eschericia  coli.  His-cdc34  was  purified  by  Ni-NTA 
chromatography  followed  by  ion  exchange  and  size-exclusion 
chromatography.  The  purity  and  concentration  of  His-cdc34  were 
determined  by  SDS-PAGE  analysis. 

In  Vitro  Ubiquitination  of  p27  and  Trfl 

Mouse  p27,  cloned  into  pCS2,  was  translated  in  vitro  in  a 
reticulocyte  lysate  system  (Promega,  Madison,  WI)  in  the  presence 
of  [35S]  -labeled  methionine.  p27  was  phosphorylated  by  purified 
recombinant  Cdk2-CyclinE  as  outlined  by  Ungermannova  et  al 
[33].  5  pi  of  the  phosphorylation  reaction  was  incubated  with  a 
ubiquitin  mixture  containing  100  nM  UBA1,  200  nM  His-cdc34, 
100  nM  SCFskp2  E3  ligase  complex,  50  nM  Cksl,  10  pM 
ubiquitin  (Sigma  Aldrich,  St.  Louis,  MO),  10  pM  methylated 
ubiquitin  (Boston  Biochem,  Cambridge,  MA),  1  pi  of  energy 
regeneration  system  (noted  as  20xER  and  consisting  of  10  mM 
ATP,  20  mM  Tris-HCl  pH  7.4,  100  mM  MgCl2,  200  mM 
creatine  phosphate,  2  mg/ml  creatine  phosphokinase,  and  10% 
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glycerol)  ,  1  |0,M  ubiquitin  aldehyde,  and  1  |tM  MG  132  in  a  total 
volume  of  1 5  (xl.  The  reaction  was  quenched  after  30  minutes  in  a 
30°C  water  bath  by  addition  of  4x  SDS  sample  bulfer.  The 
products  of  ubiquitination  were  resolved  by  SDS-PAGE,  destained 
in  a  45%  methanol  and  10%  acetic  acid  solution  in  water,  dried 
and  exposed  overnight  to  a  phosphoimager  screen,  and  scanned 
using  a  Typhoon  scanner  9400  (GE  Healthcare,  Piscataway,  NJ). 
In  vitro  Trfl  ubiquitination  was  carried  out  as  described  in  Zeng  et 
al  [34],  Recombinant  Trfl,  labeled  with  [y-33P]-ATP  by  CDK1- 
CyclinB,  was  incubated  with  0.5  pM  UBA1,  5  pM  UbcH5a, 
1  pM  SCF11’*4  E3  ligase  complex,  5  pM  ubiquitin,  100  pM 
methylated  ubiquitin,  1  pM  ubiquitin  aldehyde,  and  1  pi  20xER 
for  two  hours  at  30°C.  Ubiquitinated  Trfl  was  analyzed  by  SDS- 
PAGE  followed  by  autoradiography. 

E1/E2  Thioester  Bond  Formation  Assay 

40  nM  ubiquitin  El  (UBA1)  or  1  pM  S.  pombe  El  (Ubalp,  gift 
from  Chris  Lima)  or  0.5  pM  human  SUMO  El  (Boston  Biochem, 
Cambridge,  MA)  were  pre-incubated  with  20xER  at  30°C  for 
5  minutes  in  thioester  reaction  buffer  (20  mM  Tris  pH  7.6, 
50  mM  NaCl,  and  10  mM  MgCl2).  After  5  minutes,  1  pM 
fluorescein-ubiquitin  (Boston  Biochem,  U-590)  was  added  to 
initiate  attachment  of  ubiquitin.  All  components  were  allowed  to 
react  for  another  5  minutes  in  a  total  volume  of  5  pi.  The  reaction 
was  stopped  with  10  pi  of  SDS-PAGE  loading  buffer,  minus  DTT, 
and  the  proteins  were  resolved  using  12%  gels  that  were  run  on  ice 
to  prevent  the  reduction  of  E 1 -ubiquitin  due  to  the  heat  generated 
during  electrophoresis.  Thioester  bond  formation  was  visualized 
by  scanning  the  gel  using  Typhoon  scanner  9400  (GE  Healthcare) 
that  was  set  to  fluorescence  mode  (532  nM).  When  necessary 
100  nM  of  E2  (Cdc34)  was  added  after  the  El  enzyme  was  pre¬ 
charged  with  ATP.  Serially  diluted  Largazole  and  its  analogs  were 
incubated  with  the  reagents  as  stated  in  the  text.  ImageJ  was 
utilized  to  quantify  the  fluorescence  signal,  and  the  dose  response 
curves  were  generated  by  nonlinear  least  regression  analysis  of 
data  using  Prism  (GraphPad,  San  Diego,  CA). 

[oc-32P]-AMP:  [a-32P]-ATP  and  [32P]-PPi:[y-32P]-ATP 
Exchange  Assays 

The  reaction  mixture  contained,  in  a  final  volume  of  10  pi, 
50  mM  Tris-HCl  pH  7.5,  150  mM  NaCl,  10  mM  MgCl2 
(reaction  assay  buffer),  150  nM  human  ubiquitin  El  (UBA1), 
100  pM  ATP,  2  mM  AMP,  1  pM  [a-32P]-  or  [y-32P]-  ATP 
(Perkin  Elmer,  Waltham,  MA),  500  pM  PPi  (sodium  salt).  A  total 
of  5  pM  ubiquitin  was  added  to  the  mixture  to  initiate  the 
ATP:AMP  exchange.  After  incubation  at  30°C  for  10  minutes,  the 
reactions  were  quenched  with  EDTA,  and  0.5  pi  aliquots  of  the 
reaction  mixtures  were  spotted  on  Baker-flex®  thin  layer 
chromatography  (TLC)  polyethylenimine-modified  cellulose  plates 
(J.T.  Baker,  Phillipsburg,  NJ)  and  developed  in  filtered  0.34  M 
potassium  phosphate  pH  7.0  for  20  minutes  in  a  glass  jar.  The 
TLC  plates  were  allowed  to  air  dry  for  10  minutes,  covered  in 
plastic  wrap,  and  then  exposed  to  a  phosphoimager  plate  for  5— 
10  minutes.  The  separation  of  radiolabeled  nucleotides  was 
visualized  using  a  Typhoon  scanner  9400  (GE  Healthcare, 
Piscataway,  NJ). 

Results 

Largazole  stabilizes  GFP-p27  expression  in  Kipl 6  cells 

A  hallmark  of  many  advanced  cancers  is  an  excessive 
degradation  of  the  cyclin-dependent  kinase  inhibitor  p27,  which 
is  chiefly  directed  by  SCFskp2-mediated  ubiquitination.  Hence, 
stabilization  of  p27  degradation  represents  a  rational  approach  in 
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cancer  therapeutics.  To  identify  small  molecule  inhibitors  that  can 
stabilize  p27Kipl,  we  performed  a  screen  of  ~3000  compounds 
front  NCI  DTP  diversity  set  along  with  several  natural  products  in 
our  collection.  For  the  cell-based  screen,  we  generated  a  mink  lung 
epithelial  cell  line  (Kip  16)  stably  expressing  p27  that  was  cloned  in 
frame  with  green  fluorescent  protein  (GFP).  The  resulting  N- 
terminal  GFP-p27  fusion,  detectable  by  fluorescence  microscopy, 
was  used  to  determine  the  levels  of  p27  expression  upon  treatment 
of  cells  with  the  compound  libraries  in  96-well  format.  Much  to 
our  surprise,  the  most  potent  hit  that  emerged  from  this  screen  was 
the  natural  compound  Largazole  (L)  (Figure  1),  which  was  first 
described  by  Luescli  and  coworkers  [21]  and  subsequently 
synthesized  in  several  laboratories  including  ours  [2 1 ,24 — 
26,29,30,32] .  Largazole  induced  a  robust  and  highly  uniform 
upregulation  of  GFP-p27  at  concentrations  as  low  as  1  nM 
(Figure  2A).  As  expected,  treatment  with  the  proteasome  inhibitor 
MG  132  is  highly  effective  in  prevention  of  p27  degradation.  We 
did  not  observe  an  increase  in  GFP-p27  expression  upon 
treatment  with  the  vehicle  DMSO.  This  result  suggests  that 
Largazole  can  stabilize  GFP-p27  expression  in  cultured  cells. 

Largazole  and  select  analogs  inhibit  the  in  vitro 
ubiquitination  of  p27  and  Trfl 

Before  Largazole’s  function  as  an  inhibitor  of  histone  deacety- 
lase  was  revealed,  our  initial  investigation  into  the  mechanism  of 
this  compound  showed  its  ability  to  impede  degradation  of  GFP- 
p27  in  Kip  16  cells.  One  way  to  stabilize  p27  is  to  block  its 
ubiquitination.  Hence  we  hypothesized  that  Largazole  stabilizes 
p27  by  inhibiting  the  ubiquitination  pathway  [7,8].  One  of  the 
downsides  of  cell-based  assays  is  that  the  effects  observed  may  be 
attributed  to  the  influence  of  multiple  pathways.  For  example, 
inhibiting  the  proteasome,  elevating  transcription  of  p27,  or 
inhibiting  Cdk  activity  can  also  lead  to  an  increase  in  p27 
expression.  To  tease  out  the  mechanism  and  action  of  Largazole 
on  p27  stabilization,  we  decided  to  test  the  effect  of  Largazole  on 
p27  ubiquitination  in  a  fully  reconstituted  system  in  vitro  [33,35]. 
To  test  if  Largazole  affects  p27  ubiquitination  in  vitro,  we  added 
Largazole  to  a  p27  ubiquitin  ligation  reaction.  As  shown  in 
Figure  2C,  adding  Largazole  significantly  reduced  polyubiquiti- 
nated  p27,  suggesting  that  Largazole  blocks  p27  ubiquitination. 
Since  Largazole  is  known  to  be  a  histone  deacetylase  inhibitor  and 
has  a  thioester  moiety  that  links  an  aliphatic  chain  to  the  core,  we 
decided  to  test  whether  inhibition  of  p27  degradation  can  be 
linked  to  its  histone  deacetylase  inhibitory  activity.  The  structure- 
activity  relationship  for  Largazole  is  relatively  well  understood 
[36].  Therefore  we  next  tested  a  series  of  Largazole  analogs 
(Figure  1)  to  study  the  effect  of  structure-activity  relationship  on 
p27  ubiquitination.  To  investigate  this,  Largazole  ester  (E), 
Largazole  ketone  (K),  Largazole  macrocycle  (M),  and  seco- 
Largazole  (S)  were  tested  in  an  in  vitro  p27  ubiquitination  assay 
(Figure  2C).  We  also  added  the  HDAC  inhibitor  Trichostatin  A 
(TSA),  the  structure  of  which  can  be  found  in  Figure  1 ,  to  the  assay 
to  determine  whether  or  not  other  HDAC  inhibitors  affect  p27 
ubiquitination.  We  observed  that  Largazole  (L),  Largazole  ketone 
(K),  and  Largazole  ester  (E)  inhibited  the  ligation  of  ubiquitin  onto 
p27;  however,  the  M  and  S  analogs  and  TSA  failed  to  inhibit  the 
ubiquitination  of  p27  (Figure  2C).  The  fact  that  M  had  no 
inhibitory  activity  highlights  the  role  of  the  octanoyl  chain  in 
hindering  p27  polyubiquitination.  Seco- Largazole  (S)  did  not  affect 
p27  ubiquitination,  indicating  the  importance  of  the  topology  of 
the  inhibitory  compound.  Furthermore,  the  result  also  suggests 
that  the  thioester  moiety  of  Largazole  is  not  required  for 
inhibition,  because  the  ketone  and  ester  analogs  were  equally 
potent  in  blocking  p2 7  ubiquitination.  In  addition,  El  inhibition  is 
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Figure  1.  Chemical  structures  of  Largazole,  synthetic  analogs,  and  Trichostatin  A.  Largazole  (L)  includes  a  substituted  4-methythiazoline 
linearly  fused  to  a  thiazole,  a  3-hydroxy-7-mercaptohept-4-enoic  acid,  a  thioester  moiety,  and  a  hydrocarbon  tail.  Analogs  include  a  substituted 
ketone  (K)  and  ester  (E)  in  place  of  the  thioester  moiety,  a  macrocycle  lacking  the  thioester  moiety  and  hydrocarbon  tail  (M),  an  analog  containing  a 
macrocycle  broken  at  carbon-3  of  the  enoic  acid  (S),  and  a  thiol  analog  lacking  the  thioester  moiety  (T).  Trichostatin  A  (TSA)  contains  a  hydroxamic 
acid  functional  group. 
doi:1 0.1 371  /journal. pone.0029208.g001 


unrelated  to  HDAC  inhibitor  activity  of  Largazole  as  both  ketone 
and  ester  fail  to  inhibit  HDAC.  Prior  to  ubiquitination,  p27  is 
phosphorylated  by  the  Cdk2-CyclinE  complex.  We  carried  out  an 
in  vitro  p27  phosphorylation  assay  (as  described  in  [33])  in  the 
presence  of  either  DMSO  or  Largazole  in  order  to  test  whether  or 
not  the  decrease  in  p27  ubiquitination  was  due  to  the  inhibition  of 
the  Cdk2-CyclinE  complex.  We  observed  that  Largazole  does  not 
inhibit  the  phosphorylation  of  p27  (Figure  2B);  therefore,  the 
inhibition  of  p2  7 -ubiquitin  conjugation  is  due  to  an  inhibition  of 
the  ubiquitination  process  rather  than  phosphorylation  step. 

To  study  the  specificity  of  largazole’s  inhibition,  we  set  up  Trfl 
in  vitro  polyubiquitination  in  the  presence  of  varying  concentrations 
of  largazole  ester  and  found  that  E  inhibited  the  ubiquitin 
attachment  in  a  dose-dependent  manner  (Figure  2D).  Since 
ubiquitination  of  both  proteins  was  impeded,  and  given  that  both 
reactions  require  different  factors  to  execute  it  (p27  has  to  be 
phosphorylated  by  CDK2-CyclinE  while  there  is  no  requirement 
for  Trfl  phosphate  addition,  E2  for  p27  is  Cdc34,  while  Trfl 
needs  Ubc5Ha,  SCFSkp2  ligase  works  with  p27  while  Fbx4  is  the 
substrate  recognition  subunit  for  Trfl)  it  was  evocative  that 
largazole  compounds  stall  a  step  that  is  common  to  both 
polyubiquitination  reactions. 


Largazole  and  ester/keto  analogs  inhibit  ubiquitin  El 
activation 

Since  both  p27  and  Trfl  can  be  ubiquitinated  in  the  presence 
of  UBA1,  we  hypothesized  that  the  inhibitory  activity  of 
Largazole  is  due  to  the  deactivation  of  El.  To  test  this  hypothesis, 
we  incubated  Largazole  and  Largazole  ester  with  recombinant 
El  prior  to  carrying  out  an  in  vitro  thioester  assay  we  described 
previously  [37].  In  addition,  we  tested  the  active  thiol  form  of 
Largazole  (T)  for  El  inhibition.  The  presence  of  a  fluorescence 
signal  in  the  thioester  assay  suggests  the  formation  of  E 1  -ubiquitin 
adducts.  The  dose  dependent  decrease  in  fluorescence  indicates 
that  Largazole  and  Largazole  ester  inhibit  the  formation  of  El- 
ubiquitin  adducts  (Fig.  3 AC).  The  dose-response  curves  generated 
(data  not  shown)  suggest  an  IC50  of  approximately  29  jlM  and 
25  |iM  for  Largazole  and  Largazole  ester,  respectively.  Interest¬ 
ingly,  the  active  thiol  form  of  Largazole  (T)  failed  to  inhibit  E 1 
(Fig.  3G),  suggesting  again  that  the  octanoyl  residue  is  important 
for  inhibition. 

Activated  ubiquitin  is  normally  transferred  to  ubiquitin 
conjugating  enzymes  (E2).  If  El  activity  is  inhibited,  we  expect 
to  see  that  defects  in  E 1  activation  should  impair  the  attachment  of 
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Figure  2.  Largazole  stabilizes  p27  expression  in  Kip16  cells  and  inhibits  p27  ubiquitination  in  vitro.  (A)  Fluorescent  and  corresponding 
bright-field  images  of  Kipl  6  cells  treated  with  varying  concentrations  of  Largazole  (L).  L  treatment  induces  the  expression  of  GFP-p27  in  a  dose- 
dependent  fashion.  Addition  of  MG1 32  (1  pM)  prevents  the  degradation  of  GFP-p27  via  the  ubiquitination  and  subsequent  proteasomal  degradation 
pathway.  The  vehicle  control,  DMSO,  has  no  effect  on  the  reporter  protein  stabilization.  (B)  L  fails  to  inhibit  the  phosphorylation  of  p27  by  the  Cdk2/ 
CyclinE  complex  compared  to  the  positive  control.  L  (250  pM,  lane  3,  and  125  pM,  lane  4)  was  incubated  with  the  Cdk2/CyclinE  complex  prior  to  the 
autophosphorylation  of  Cdk2/CyclinE  step.  Phosphorylated-p27  was  identified  by  protein  standard.  (C)  L,  K,  and  E  reduce  polyubiquitinated  forms  of 
p27  while  M  and  S  have  no  inhibitory  effects.  Ubiquitin-activating  enzyme  El  (100  nM),  UBA1,  was  incubated  with  100  pM  of  each  compound  prior  to 
the  reaction.  (D)  E  reduces  polyubiquitinated  forms  of  T rf  1  in  a  dose-dependent  fashion.  UBA1  (100  nM)  was  incubated  with  either  DMSO  or  various 
concentrations  of  E  ranging  from  250  pM  to  1  pM  prior  to  the  reaction. 
doi:10.1371/journal.pone.0029208.g002 


ubiquitin  onto  Cdc34  (E2).  To  further  validate  El  inhibition,  we 
included  Cdc34,  the  E2  enzyme  required  for  p27  ubiquitinaton,  in 
the  El  reaction  mixture.  As  shown  in  Figure  3EF,  in  the  presence 
of  ATP,  fluorescent  ubiquitin  is  transferred  to  Cdc34  indicated  by 
the  presence  of  a  fluorescent  Cdc34  band  on  the  gel.  Upon 
incubation  with  El,  Largazole  or  Largazole  ester  reduce  the 
amount  of  ubiquitin  molecules  that  are  transferred  from  El  to  E2 
in  a  dose-dependent  fashion  (Fig.  3BD). 

The  decreased  ubiquitin  transfer  could  be  attributed  to  either 
El  or  E2  inhibition;  therefore,  we  produced  El  precharged  with 
ubiquitin  by  incubating  ATP  and  fluorescent  ubiquitin  for 
15  minutes  at  room  temperature  followed  by  the  addition  of 
Cdc34,  which  was  preincubated  with  either  Largazole  or 
Largazole  ester.  If  either  compound  inhibits  the  transfer  of 
ubiquitin  from  El  to  E2,  then  we  would  observe  a  significant 
decrease  in  Cdc34  fluorescence  regardless  of  the  order  we  added 
the  compounds.  Interestingly,  Largazole,  preincubated  with 
Cdc34,  fails  to  inhibit  the  transfer  of  ubiquitin  from  precharged 
El  at  concentrations  <1  mM  (Fig.  3E).  Furthermore,  in  a  similar 
experiment,  Largazole  ester  begins  to  inhibit  the  transfer  of 
ubiquitin  from  precharged  El  to  Cdc34  at  concentrations  around 
500  pM  (Fig.  3F),  although  this  concentration  is  significantly 
above  the  IC5o  of  El  inhibition.  These  results  suggest  that 
Largazole  and  Largazole  ester  exhibit  selectivity  towards  ubiquitin 
E 1 .  Also,  this  result  suggests  that  either  compound  fails  to  promote 
the  hydrolysis  of  ubiquitin  thioesters  on  precharged  El. 


Largazole  ketone  inhibits  the  adenylation  step  of  El 
activation 

El  forms  an  ubiquitin-adenylate  intermediate  during  the 
course  of  its  catalytic  cycle  [3] .  Thus  the  mechanism  of  ubiquitin 
El  activation  can  be  studied  by  assaying  ATP:PPi  and  ATP:AMP 
exchanges  [3].  Production  of  AMP  in  the  [a-3‘P]-AMP:[  a-32P]- 
ATP  exchange  assay  guarantees  that  a  thioester  bond  is  formed 
between  El  and  ubiquitin,  while  the  release  of  PPi,  measured  by 
the  [32P]-PPi:[y-32P]-ATP  exchange  assay,  signals  the  formation 
of  ubiquitin  adenylate.  To  further  dissect  the  mechanism  of 
Largazole  inhibition,  two  nucleotide  exchange  assays  were 
carried  out  in  the  presence  of  Largazole  derivatives.  For  these 
experiments  we  used  Largazole  ketone,  which  is  similar  to 
Largazole  and  Largazole  ester.  From  the  results  shown  in 
Figure  4,  it  is  evident  that  the  first  two  concentrations  of 
Largazole  ketone  (100  and  50  pM)  inhibit  ubiquitination  of  El 
similarly  and  were  also  inhibitory  in  both  types  of  exchange 
assays.  The  lack  of  a  [3“P]-PPi  signal  suggests  that  the  adenylation 
step  did  not  occur;  consequently,  ubiquitin  could  not  be 
transferred  to  the  active  site  cysteine  to  trigger  the  release  of 
AMP.  Both  steps  of  the  El -catalyzed  reactions  can  be  measured 
by  the  AMP:ATP  exchange  assay.  The  lack  of  an  [a-3“P]-AMP 
signal  further  suggests  that  the  adenylation  step  is  inhibited  by 
Largazole  ketone.  Thus  Largazole  or  Largazole  derivatives  act  on 
the  first  step  of  ubiquitin  activation  pathway  by  blocking  the 
formation  of  ubiquitin-adenylate. 
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Figure  3.  Largazole  (L)  and  largazole  ester  (E)  inhibit  ubiquitin  El  in  a  dose  dependent  manner  in  vitro.  (A,C)  L  and  E  inhibit  transfer  of 
ubiquitin  onto  El  in  a  concentration-dependent  manner.  Thioester  assay  of  El  activity  using  fluorescein  ubiquitin  (Ub-F).  Thioester  bond  formation 
between  El  and  Ub-F  is  ATP-dependent  (lane  2  vs.  lane  1).  In  addition,  DMSO  has  no  effect  on  the  formation  of  the  thioester  linkage  as  seen  in  lane  2 
of  both  gels.  50  nM  El  was  incubated  with  decreasing  concentrations  of  L  (A)  or  E  (C)  for  15  minutes  at  room  temperature  followed  by  addition  of  a 
cocktail  containing  ATP  and  Ub-F.  After  5  minutes  of  incubation,  the  reactions  were  resolved  by  SDS-PAGE  under  non-reducing  conditions.  Ub-F  was 
used  to  show  equal  loading.  (B,D)  Thioester  assay  of  the  ubiquitin  transfer  from  El  to  E2  (Cdc34).  Largazole  or  Largazole  ester,  when  preincubated 
with  50  nM  El  for  15  minutes,  inhibit  the  transfer  of  ubiquitin  from  El  to  Cdc34  in  a  concentration-dependent  manner.  (E)  Largazole  selectively 
inhibits  the  activity  of  El  not  E2.  50  nM  El  was  pre-charged  with  ATP  and  then  added  to  Cdc34  that  was  previously  incubated  with  decreasing 
concentrations  (1  mM-16  pM)  of  L  in  thioester  reaction  mixture.  (F)  Largazole  ester  inhibits  E2  at  high  concentrations.  Pre-charged  El  was  added  to 
reactions  that  contained  Cdc34  pre-incubated  with  E  ranging  from  1  mM  to  16  pM  and  resolved  by  SDS-PAGE  under  non-reducing  conditions. 
Complete  inhibition  of  ubiquitin  transfer  to  E2  was  observed  at  1  mM  of  E,  with  only  modest  inhibition  at  500  pM.  (G)  Largazole  thiol  (T)  has  no  effect 
on  transfer  of  ubiquitin  onto  El.  The  reaction  was  carried  out  as  described  in  A,C. 
doi:1 0.1 371  /journal.pone.0029208.g003 
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Selectivity  of  Largazole  ketone  against  SUMO  El  and 
Ubalp 

In  addition  to  ubiquitin,  there  exist  several  ubiquitin-like 
proteins  that  covalently  modify  other  proteins.  All  of  die 
ubiquitin-like  proteins  have  activation  pathways  similar  to 
ubiquitin  [38].  In  order  to  study  the  specificity  of  Largazole  to 
the  ubiquitin  pathway,  we  incubated  Largazole  ketone  with 
SUMO-activating  El  enzyme  prior  to  carrying  out  a  thioester 
assay.  From  the  results  in  Figure  5B,  we  found  that  Fargazole 
ketone  is  ineffective  in  inhibiting  the  formation  of  El -SUMO 
adducts.  From  the  dose-response  curve  generated  from  the  SUMO 
El  fluorescence  results,  the  IC50  is  approximately  450  |U,M  as 
opposed  to  30  pM  for  ubiquitin  El  (data  not  shown).  Thus 
Largazole  is  relatively  selective  in  perturbing  ubiquitin  El 
activation. 

Ubiquitin  and  the  ubiquitin  El  enzyme  are  highly  conserved 
among  eukaryotes  [38].  Sequence  analysis  shows  a  45%  homology 
between  the  human  ubiquitin-activating  enzyme  El  (UBA1)  and 
S.pombe  El  (ptr3 /Ubalp)  at  the  amino  acid  sequence  level.  To  test 
whether  Largazole  ketone  inhibits  the  S.  pombe  E 1 ,  we  carried  out  a 
thioester  assay  using  Largazole  ketone  and  die  ubiquitin  El 


homologue  in  S.  pombe,  Ubalp.  The  results  in  Figure  5A  suggest 
that  Largazole  ketone  fails  to  inhibit  the  formation  of  E 1  -ubiquitin 
adducts  at  concentrations  less  than  1  mM.  Taken  together,  these 
results  suggest  that  Largazole  and  its  derivative  are  highly  selective 
in  inhibiting  the  ubiquitin  El  enzyme. 

Discussion 

In  this  study,  we  showed  that  Largazole  and  its  analogs 
selectively  inhibit  ubiquitin  El  enzyme  activity  in  vitro.  Also,  we 
demonstrated  that  the  inhibitory  activity  of  Largazole  is 
independent  of  its  inhibitory  activity  towards  the  histone 
deacetylase  enzymes.  Structure-activity  relationship  analysis  shows 
that  the  thioester  bond  is  not  required  for  inhibition  but  the 
macrocycle  core  and  aliphatic  tail  are  indispensible.  Largazole 
blocks  ubiquitin  activation  at  the  adenylation  step  and  without 
perturbing  ubiquitin  transfer  from  El  to  E2.  Finally  we  show  that 
Largazole  inhibition  of  E 1  is  highly  selective  as  it  does  not  inhibit  a 
highly  related  ubiquitin  El  enzyme  from  S.  pombe  and  is  almost 
twenty  fold  less  effective  in  inhibiting  die  activation  of  SUMO  E 1 . 
Taken  together,  our  results  reveal  that  Largazole  represents  a  new 
class  of  ubiquitin  El  inhibitors. 
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Figure  4.  Largazole  ketone  inhibits  the  adenylation  of  the  El 
ubiquitin-activating  enzyme.  (A)  Largazole  ketone  inhibits  ligation 
of  ubiquitin  onto  El  in  a  concentration-dependent  fashion.  Reduction 
of  El  ~Ub  adducts  was  determined  by  thioester  assay  utilizing 
fluorescein  ubiquitin.  (B)  Largazole  ketone  inhibits  the  adenylation 
step  in  ubiquitin  El  activation  in  a  concentration-dependent  fashion.  K 
was  serially  diluted  (100  pM  to  1.5  pM)  and  incubated  with  UBA1 
(150  nM)  at  room  temperature  for  five  minutes.  The  thioester  reaction 
mixture  was  mixed  with  ubiquitin  to  initiate  the  PPi:ATP  exchange 
(middle  panel)  or  AMP:ATP  exchange  (bottom  panel)  and  added  to  the 
UBA1/K  mixture.  All  reactions  were  halted  with  addition  of  EDTA  after 
10  minute  incubation  at  37  C,  resolved  using  Cellulose  PEI  TLC  plates, 
and  analyzed  using  a  phosphoimager. 
doi:1 0.1 371  /journal.pone.0029208.g004 


We  identified  that  Largazole  caused  a  robust  increase  in  GFP- 
p27  expression  in  Kip  16  cells.  This  observation  led  us  to  further 
investigate  the  mechanism  of  GFP-p27  stabilization  by  Largazole. 
Using  an  in  vitro  ubiquitination  assay,  we  were  able  to  delineate  the 
inhibitory  point  where  Largazole  acts  on  in  the  ubiquitination 
pathway,  namely  the  El  enzyme.  However,  there  is  a  disconnect 
between  the  potency  of  El  inhibition  in  vitro  and  GFP-p27 
stabilization  in  cells.  The  EC50  of  Largazole  for  GFP-p27 
stabilization  is  in  the  low  nM  range,  yet  El  inhibition  is  at 
~30  pM.  This  results  suggests  that  the  stabilization  of  GFP-p27  is 
unlikely  caused  by  El  inhibition,  but  is  most  likely  a  result  of 
HDAC  activity,  which  is  known  to  block  cell  cycle  progression  and 
cause  cell  growth  arrest.  Consistent  with  this  hypothesis,  Largazole 
ketone  and  ester,  two  Largazole  analogs  that  do  not  inhibit 
HDACs,  do  not  increase  GFP-p27  levels  when  Kip  16  cells  were 
treated  (data  not  shown).  However,  other  interpretations  may 
account  for  the  failure  of  Largazole  ketone  or  ester  to  raise  GFP- 
p27  by  inhibiting  El  in  cells.  For  example,  we  do  not  know  if  or 
how  these  compounds  penetrate  cells  and  how  stable  they  are  once 
they  enter  the  cells.  These  investigations  have  to  be  undertaken 
before  these  analogs  can  be  further  developed  for  in  vivo 
applications. 


Figure  5.  Investigation  into  the  selectivity  of  Largazole  ketone. 

(A)  Largazole  ketone  (K)  fails  to  inhibit  the  ligation  of  ubiquitin  onto 
Ubalp,  a  homologue  of  UBA1  from  S.  pombe.  Formation  of  Ubalp- 
ubiquitin  adducts  was  determined  by  thioester  assay  utilizing 
fluorescein-ubiquitin.  Ubalp  (1.03  pM)  was  incubated  with  either 
DMSO  or  various  concentrations  of  K  serially  diluted  from  1000  pM  to 
31  pM.  (B)  K  inhibits  ligation  of  SUMO-1  onto  human  SUMO  El  in  a 
concentration-dependent  fashion.  Reduction  of  E1-SUMO  adducts  was 
determined  by  thioester  assay  utilizing  fluorescein-SUMO-1.  hSUMO  El 
(500  nM)  was  incubated  with  either  DMSO  or  various  concentrations  of 
K  serially  diluted  from  1000  pM  to  31  pM. 
doi:1 0.1 371/journal.pone.0029208.g005 

Panepophenanthrin,  a  natural  compound  derived  from  the 
mushroom  strain  Panus  rudis,  and  Himeic  acid  A,  derived  from  the 
marine  fungus  Aspergillus,  are  the  first  and  second  discovered 
inhibitors  of  the  ubiquitin-activating  enzyme  El,  respectively 
[39,40],  Both  compounds  were  tested  in  vitro  using  recombinant 
El;  however,  the  cellular  activity  and  mechanism  were  not 
determined  [41].  PYR-41  and  related  pyrazones  are  another  set  of 
compounds  that  were  discovered  to  inhibit  ubiquitin  El  and  the 
first  set  of  E 1  inhibitors  described  to  enter  cells  and  differentially 
kill  transformed  cells  [42].  The  IC5q  of  PYR-41  is  around  5  pM, 
thus  more  potent  than  the  compounds  described  here.  However, 
the  exact  mechanism  of  PYR-41  inhibition  is  not  known.  Ub- 
AMSN  represents  a  distinct  class  of  protein  based  inhibitors  of 
ubiquitin  El.  Ub-AMSN  contains  a  sulfamide  group  attached  to 
the  carboxyl  terminus  of  ubiquitin  as  a  nonhydrolyzable  mimic  of 
the  phosphate  group  in  the  cognate  Ub/Ubl-AMP  adenylate 
intermediate.  Thus,  like  Largazole  analogs,  it  blocks  the  first  step 
of  El  reaction  [43,44],  Unfortunately,  Ub-AMSN  cannot  be  used 
in  cells  as  it  cannot  pass  through  the  cell  membrane.  However, 
Ub-AMSN  turns  out  be  a  very  useful  for  probing  the  structure  and 
biochemical  mechanisms  of  El  enzyme  [44],  Therefore,  Largazole 
and  analogs  could  also  be  useful  tools  for  probing  ubiquitin 
function. 

One  of  the  most  important  questions  to  be  answered  is  whether 
or  not  ubiquitin  or  ubiquitin-like  El  inhibitors  are  therapeutically 
relevant.  Since  only  one  ubiquitin  El  enzyme  is  responsible  for  a 
majority  of  protein  ubiquitination  in  humans,  inhibiting  El  will 
influence  the  degradation  of  proteins  across  several  pathways  and 
may  lead  to  toxicity  and,  consequently,  poor  therapeutic  efficacy. 
Bortezomib  is  die  first  FDA-approved  proteasome  inhibitor  for  the 
treatment  of  relapsed/ refractory  myeloma  and  mantle  cell 
lymphoma.  [19].  The  proteasome,  particularly  the  26S  protea¬ 
some,  is  the  final  step  in  ubiquitin-mediated  protein  degradation 
and  regulates  various  pathways  necessary  for  cellular  function. 
The  clinical  success  and  efficacy  of  Bortezomib  gives  rise  to  the 
possibility  that  inhibitors  of  ubiquitin  El  will  also  share  similar 
success.  NEDD8  is  a  protein  modifier  that  shares  mechanistic  and 
structural  similarities  to  ubiquitin.  Currently,  the  cullin  family  of 
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proteins  has  been  characterized  as  the  target  for  NEDD8 
conjugation  [11].  MLN4924  is  a  potent  and  selective  inhibitor 
of  the  NEDD8-activating  enzyme  (NAE)  that  exhibited  potent 
cytotoxicity  against  several  human  tumor-derived  cell  lines  [45]. 
Interestingly,  MLN4924  shares  a  similar  mechanism  to  Largazole 
analogs.  MLN4924  reacts  covalently  with  NEDD8  mimicking  a 
NEDD8  adenylate  that  is  incapable  of  driving  the  reaction 
forward,  therefore,  blocking  the  activity  of  NAE  [46].  MLN4924  is 
currently  undergoing  phase  I  clinical  trials  in  patients  with 
lymphoma,  multiple  myeloma,  or  any  form  of  nonhematologic 
malignancies.  The  in  vitro  and  possible  clinical  success  of  the  NAE 
inhibitor  MLN4924  further  supports  the  concept  that  El 
inhibitors  are  potential  promising  cancer  therapeutics. 

Our  preliminary  structure  activity  relationship  studies  suggest 
that  the  pro-drug  form  of  Largazole  including  both  the 
hydrocarbon  tail  and  the  macrocycle  are  essential  for  El 
inhibition.  For  Largazole  analogs  to  be  developed  as  potential 
antitumor  drugs,  additional  analogs  are  needed  to  be  synthesized 
in  order  to  improve  its  potency  toward  ubiquitin  El.  The  most 
promising  aspect  of  Largazole  analogs  as  ubiquitin  E 1  inhibitors  is 
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the  selectivity  and  specificity  of  Largazole.  Largazole  analogs  not 
only  display  discrimination  over  related  SUMO  El  enzyme  but 
also  remarkable  selectivity  in  targeting  human  ubiquitin  El. 
Future  structural  studies  would  be  helpful  to  understand  how 
Largazole  analogs  inhibit  El,  and  insights  gained  front  such 
studies  may  help  to  develop  more  specific  inhibitors  of  El. 
Experiments  to  test  these  hypotheses  are  currently  underway. 
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17  Abstract 

18  Understanding  how  cells  migrate  individually  and  collectively  during  development  and 

19  cancer  metastasis  can  be  significantly  aided  by  a  computation  tool  to  accurately  measure  not 

20  only  cellular  migration  speed,  but  also  migration  direction  and  changes  in  migration  direction  in 

21  a  temporal  and  spatial  manner.  We  have  developed  such  a  tool  for  cell  migration  researchers, 

22  named  Pathfinder,  which  is  capable  of  simultaneously  measuring  the  migration  speed,  migration 

23  direction,  and  changes  in  migration  directions  of  thousands  of  cells  both  instantaneously  and 

24  over  long  periods  of  time  from  fluorescence  microscopy  data.  Additionally,  we  demonstrate  how 

25  the  Pathfinder  software  can  be  used  to  quantify  collective  cell  migration.  The  novel  capability  of 

26  the  Pathfinder  software  to  measure  the  changes  in  migration  direction  of  large  populations  of 

27  cells  in  a  spatiotemporal  manner  will  aid  cellular  migration  research  by  providing  a  robust 

28  method  for  determining  the  mechanisms  of  cellular  guidance  during  individual  and  collective  cell 

29  migration. 


2 


30  Introduction 

31  Cellular  migration  has  been  shown  to  be  an  important  process  in  cancer  progression, 

32  development,  tissue  repair,  and  immune  response  [1-10].  As  a  result,  a  plethora  of  research 

33  has  been  performed  to  identify  the  molecular  mechanisms  behind  how  individual  cells 

34  achieve  migration,  as  well  as  how  neighboring  cells  migrate  cooperatively  in  collective 

35  migration  (reviewed  in  [1 1  -1 3]  and  [1 4],  respectively).  Collective  migration  is  defined  as  the 

36  ability  of  physically  interacting  cells  to  adopt  a  common  migration  direction  [14,  15].  Like 

37  individual  cell  migration,  the  collective  migration  of  cells  has  been  shown  to  be  an  important 

38  process  in  cancer  progression,  development  and  wound  repair  [16-23].  Such  collective 

39  behavior  results  from  each  cell  responding  to  the  environmental  stimuli  of  neighboring  cells, 

40  in  addition  to  non-cell  environmental  stimuli  [4,  5,  14,  15,  17,  19,  20,  24-32].  Although  a 

41  relatively  large  amount  of  research  has  been  conducted  to  determine  mechanisms  behind 

42  individual  cell  migration,  far  less  is  known  about  exactly  how  cells  migrate  collectively. 

43  Furthermore,  there  is  no  standard  method  in  the  literature  to  quantify  the  ‘collectiveness’ 

44  behavior  during  collective  migration  [33-35]. 

45  Previous  research  into  individual  cell  migration  has  revealed  important  fundamental 

46  mechanisms  by  which  cells  migrate.  For  instance,  when  an  individual  cell  migrates  on  a  two- 

47  dimensional  (2D)  surface,  it  projects  a  front  end  extension  that  can  either  be  broad  (termed 

48  a  llamelipodia)  or  with  multiple  spike-like  extensions  (termed  fillipodia),  which  are  the  result 

49  of  coordinated  polymerization,  depolymerization,  and  branching  of  the  actin  cytoskeleton 

50  [12,  24,  36-46].  Such  coordination  of  actin  dynamics  is  controlled  by  local  recruitment  of  cell 

51  polarity  maintain  proteins,  such  as  CDC42/Rac  and  Rho,  which  either  directly  or  indirectly 

52  regulate  actin  structure,  polymerization,  and  attachment  to  the  extracellular  matrix  [37,  38, 

53  47-51].  The  attachment  of  the  actin  cytoskeleton  is  largely  mediated  by  protein  complexes, 

54  termed  focal  adhesions,  which  anchor  the  actin  cytoskeleton  to  trans  membrane  integrin 
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55  receptors  and  the  extracellular  matrix  [52],  The  assembly  of  focal  adhesions  allows  for  the 

56  cell  to  successfully  attach  a  front  end  extension  to  the  extracellular  matrix  and  the 

57  disassembly  of  focal  adhesions  allows  a  cell  to  detach  the  rear  during  rear  end  retraction 

58  [53-55].  Focal  adhesion  turnover  and  the  resulting  changes  to  the  actin  cytoskeleton  are 

59  regulated  by  several  kinase  activities,  including  focal  adhesion  kinase  (FAK),  Src  kinase  and 

60  Rho  GTPase  [12,  56-60],  The  temporal  and  spatial  regulation  of  both  the  actinomyosin 

61  skeleton  and  focal  adhesions  are  regulated  by  a  complex  combination  of  growth  factor 

62  signaling  and  extracellular  matrix  protein  activities,  which  influence  the  speed  of  actin  and 

63  focal  adhesion  dynamics,  ultimately  influencing  how  fast  a  cell  can  migrate  [61 , 62], 

64  Our  current  understanding  of  the  biochemical  mechanisms  underlying  cellular  migration 

65  have  been  primarily  the  result  of  in  vitro  studies  conducted  in  2D  cell  culture  model  systems 

66  [2,  11-13,  24,  31, 32,  38,  40-43,  54,  56,  57,  60,  62-70].  Flowever,  several  critical  biochemical 

67  activities  governing  cell  migration  have  proven  to  play  similar  roles  in  three  dimensional  (3D) 

68  model  systems  and  in  vivo.  For  instance,  in  2D,  3D  and  in  vivo  experiments,  CDC42/Rac 

69  activity  determine  cellular  polarity  [71 , 72],  Similarly,  FAK  kinase  mediates  cellular  migration 

70  both  in  2D  and  3D  assays  [73-75].  As  a  result,  investigations  performed  in  2D  assays  have 

71  shed  light  on  biochemical  mechanisms  that  have  proven  to  have  physiological  relevance. 

72  Flowever,  recent  research  has  also  revealed  that  there  is  significant  differences  in  cell 

73  migration  machinery  between  cells  in  2D  versus  3D  [75-77].  Although  the  conclusions  made 

74  in  2D  migration  studies  will  always  require  confirmation  of  physiological  relevance  in  in  vivo 

75  studies,  they  remain  a  valuable  tool  for  initial  investigations  into  the  molecular  mechanisms 

76  behind  cellular  migration  compared  to  3D  and  in  vivo  studies  because  they  allow  for  tight 

77  control  of  experimental  conditions  and  more  accurate  observation  of  cellular  migration 

78  behavior  at  single  cell  resolution  without  the  use  of  relatively  complex  microscopes,  such  as 

79  two-photon  and  confocal  microscopes.  Many  of  the  concerns  about  discrepancies  in 

80  biochemical  mechanisms  behind  2D  and  3D  motility  may  prove  to  be  overcome  by  imaging 
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81  individual  cell  motility  in  2D  on  soft  extracellular  matrices,  which  have  been  shown  to  be 

82  more  closely  similar  to  in  vivo  tissues  than  plastic  or  glass  cell  culture  plates  [78]. 

83  The  behavior  of  migrating  cells  can  be  characterized  by  migration  speed,  migration 

84  direction,  and  migration  persistence  (the  ability  of  a  cell  to  maintain  its  migration  direction). 

85  In  2D  studies,  the  measurement  of  cell  migration  behavior  is  conducted  by  either  manual 

86  cell  tracking  [79-81  ]  or  automated  cell  tracking  [31 , 32,  61 , 82-85].  Such  cell  tracking 

87  experiments  have  not  only  shed  light  on  how  a  cell  achieves  migration,  but  also  have  shown 

88  that  cells  can  undergo  chemotaxis  towards  a  localized  biochemical  signals  [20,  86,  87],  In 

89  these  studies,  a  Dunn  Chamber  is  used  to  present  a  chemokine  gradient  to  cells,  where 

90  cells  migrate  upstream  of  EGF  and  uPA  gradients  [88].  Such  studies  into  how  cells  achieve 

91  chemotaxis  highlight  the  need  for  cell  migration  tracking  programs  to  not  only  calculate  the 

92  speed  and  persistence  of  cells,  but  also  to  report  the  direction  and  changes  in  direction 

93  during  cellular  migration. 

94  Although  several  computational  tools  exist  that  allow  for  automated  cell  tracking  of 

95  individual  cells  in  time-lapse  microscopy  videos  (Fig.  1),  these  tools  focus  almost  entirely  on 

96  either  the  speed  or  persistence  (ability  to  maintain  a  migration  direction)  or  a  cell.  We  have 

97  developed  an  automated  high  throughput  cell  tracking  software,  named  Pathfinder,  which  is 

98  capable  of  simultaneously  measuring  and  reporting  cellular  migration  speed,  migration 

99  direction  and  changes  in  migration  direction  of  thousands  of  fluorescently  labeled  cells  for 

100  an  unlimited  number  of  microscopy  videos.  The  Pathfinder  software  has  two  improved 

101  features  that  distinguish  it  from  previous  cell  motility  computational  tools,  as  the  Pathfinder 

102  output  is  able  to  report  instantaneous  cell  migration  direction  and  changes  in  cellular 

103  migration  direction.  Although  the  cellular  migration  field  has  elucidated  many  fundamental 

104  mechanisms  behind  cellular  migration,  there  are  several  key  questions  that  remain  to  be 

105  answered.  Specifically,  ‘how  does  a  cell  select  a  migration  direction?’  and  ‘what  makes  a 
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106  cell  change  direction?’  Answering  these  questions  will  require  the  ability  to  study  the 

107  migration  direction  and  changes  in  migration  direction  at  single  cell  and  instantaneous 

108  resolution.  Additionally,  we  explain  how  measurements  of  cellular  migration  direction  can  be 

109  used  to  quantify  collective  migration.  In  summary,  the  Pathfinder  software  aims  to  propel  the 

110  cell  migration  field  forward  in  both  the  mechanistic  understanding  of  individual  cell  migration, 

111  as  well  as  collective  cell  migration  by  allowing  researchers  to  fully  characterize  cell  migration 

112  behavior  in  an  automated  and  high  throughput  manner. 

113 

114  Methods 

115 

116  Fluorescent  Labeling  of  Cells,  Cell  Culture  and  Cellular  Imaging 

117  Stable  transgenic  HaCaT  (Cell  Lines  Services,  Germany)  and  MDA-MB-231  (ATCC, 

118  HTB-26)  cell  lines  were  fluorescently  labeled  via  retroviral  mediated  gene  transfer  of 

119  mCherry-Histone  H2B  using  the  pRex-mCherry-H2B  plasmid.  For  all  experiments,  cells 

120  were  cultured  in  DMEM  lacking  phenolphthalein  red  and  supplemented  with  2  mM  L- 

121  glutamine,  1 00  Units/mL  penicillin  and  1 00  pg/mL  streptomycin.  For  low  density  assays, 

122  cells  were  plated  at  an  average  density  of  300  cells/mm2  for  both  PlaCaT  and  MDA-MB-231 

123  cells.  For  confluent  monolayer  experiments,  HaCaT  cells  were  plated  at  an  average  density 

124  of  1000  cells/mm2.  For  epithelial  sheet  assays,  HaCaT  cells  were  plated  at  an  average 

125  density  of  1 200  cells/mm2  for  3  hours  at  37  qC,  after  which  half  of  partially  adherent  cells 

126  were  manually  removed  using  a  200  pi  pipet  tip.  TGF(3  stimulations  were  conducted  with  100 

127  pM  ligand,  while  EGF  stimulations  were  conducted  with  100  nM  ligand.  An  ImageXpress 

128  MicroXL  high  throughput  wide-field  fluorescence  microscope  (Molecular  Devices)  was  used 

129  for  imaging  experiments  at  37  °C  and  5  %  C02.  All  microscopy  videos  were  acquired  with  a 

130  frame  rate  of  one  frame  every  seven  minutes,  in  the  mCherry  fluorescence  channel  at  a 
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magnification  of  lOx,  where  each  pixel  represents  1 .314  pm  at  a  pixel  binning  of  2  x  2.  The 
accompanying  MetaXpress  software  was  used  to  compile  video  files  from  time-lapse  images 
for  each  well  of  a  96  well  plate. 

Programming,  Input  Parameters,  and  Output  for  the  Pathfinder  Software 


The  Pathfinder  software 

(https://universitvofcolorado.box.eom/s/qzs6nos4470sifsaw8wq  and  Fig.  SI )  was 

written  in  the  Java  programming  language.  User  specification  is  required  for  cell  radius 
(pixels),  minimum  track  length,  the  interval  of  frames  for  desired  calculations  (frame  n  - 
frame  n+a,  where  a  represents  the  number  of  frames  to  skip  for  calculations),  percentage  of 
pixels  in  the  video  that  represent  cells,  and  the  directory  path  for  the  folder  containing  .avi 
files  (Fig.  S2).  The  output  for  each  video  file  is  a  single  Excel  spreadsheet  (Fig.  S3  and  S4). 
The  Pathfinder  software  requires  only  decompression  of  the  attached  .zip  file  and 
installation  of  JAVA  runtime  environment  on  either  a  32-bit  or  64-bit  Windows  Machine. 
Please  note  that  use  of  pathfinder  on  a  64-bit  machine  allows  for  higher  memory  use  in 
Java,  which  allows  for  analysis  of  greater  numbers  of  cells  in  a  single  video. 


Calculation  of  Migration  Parameters,  Persistence  Time  and  Nearest  Neighbor 
Analyses 

Cellular  speed  was  calculated  as  the  displacement  of  a  cell  (pixels)  over  1  frame. 
Conversion  to  pm/hour  is  determined  by  the  following  equation: 


Conversion  Factor  = 


( Camera  Pixel  SizexBinning  Factor ) 
C Magnification  FactorxFrame  Rate ) 


The  Angle  of  Trajectory  was  calculated  from  the  following  discontinuous  equations: 
if(dx  >  0 ,cLy>  0),  then  0Trajectory  =  360  —  ArcTAN  (^j 
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if(dx  >  0,  dy  <  0),  then 
if(dx  <  0,  dy  >  0),  then 

if(dx  <  0,  dy  <  0),  then 
if(dx  =  0  ,dy  =  0),  then 
if(dx  >  0,  dy  =  0),  then 
if(dx  <  0,  dy  =  0),  then 
if(dx  =  0,  dy  <  0),  then 
if(dx  =  0,  dy  >  0),  then 


® Trajectory  ~  —ArcTAN  ^ 
0 Trajectory  =  180  —  ArcTAN 

0 Trajectory  =  180  —  ArcTAN 
@  Trajectory  N O  Angle 

® Trajectory  0 

® Trajectory  180 
® Trajectory  90 
® Trajectory  ~  270 


Angle  of  Deflection  was  calculated  from  the  following  discontinuous  equations: 

if  (0 Deflection  —  180),  then  6 Deflection  ~  (^Trajectory  t=n+l  ~  ^Trajectory  t=n  )  —  360 
if  (0 Deflection  —  — 180 ^,then  6 Deflection  ~  (^Trajectory  t=n+l  —  ^Trajectory  t=n  )  4"  360 


Persistence  time  calculations  were  performed  using  a  modified  in-house  MatLab  based 
program  developed  by  Dr.  Douglas  Lauffenburger  (MIT)  (Fig.  S5). 

Nearest  neighbor  calculations  were  done  in  Excel  using  an  inter-centroid  distance  matrix 
of  (all  cells)  x  (all  cells)  for  each  frame.  Nearest  neighbors  were  defined  as  cells  whose 
centroids  are  within  100  pm  of  each  other.  Paired  random  migration  index  of  the  angle  of 
trajectory  (PRMI  ©Trajectory)  calculations  were  also  done  in  Excel,  where  the  standard 
deviation  of  the  migration  directions  between  pairs  of  neighboring  cells  were  averaged  over 
greater  than  100  pairs  of  cells  for  each  condition.  This  average  of  standard  deviations  is 
referred  to  as  the  PRMI  ©Trajectory 
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Results 


The  overview  and  capabilities  of  the  pathfinder  software 

The  JAVA  based  Pathfinder  software  was  developed  to  allow  researchers  to  easily 
analyze  large  data  sets  of  time-lapse  fluorescence  microscopy  videos  of  motile  cells.  Since 
cellular  tracking  is  already  a  well-established  technique,  our  software  implements  a 
previously  validated  tracking  algorithm  (‘Particle  Tracker’)  developed  by  Sbalzarini  et.  a I  to 
detect  each  fluorescently  labeled  nuclei  in  each  frame  (Fig.  2A,  left),  as  well  as  to  assemble 
such  positional  information  into  cellular  tracks  (Fig.  2A,  right),  as  described  in  their 
publication  [89].  Since  cellular  positions  alone  are  of  little  use  to  researchers  in  the  cell 
migration  field,  we  developed  an  analysis  algorithm  to  transform  the  previous  ‘Particle 
Tracker’  output  into  an  excel  spreadsheet  that  displays  calculations  of  the  speed,  the 
direction,  and  changes  in  direction  of  individual  cells,  as  well  as  the  average  values  for  a 
population  of  cells  (Fig.  2B,  S3,  and  S4).  Simultaneous  reporting  of  these  parameters 
makes  the  Pathfinder  software  unique  compared  to  other  available  computational  motility 
programs  (Fig.  1).  In  addition,  Pathfinder  is  capable  of  running  batch  parallel  processing  of 
unlimited  .avi  files,  allowing  for  automated  and  high  throughput  data  processing  of 
fluorescent  time-lapse  microscopy  videos,  provided  they  are  placed  in  a  single  folder  that 
can  be  navigated  to  from  within  the  Pathfinder  GUI.  The  Pathfinder  program  requires  either 
a  32  bit  or  a  64  bit  windows  operating  system  with  JAVA  Runtime  Environment,  and  is 
available  at  https://universitvofcolorado.box.eom/s/azs6nos4470sifsaw8wg  . 

Fluorescent  Labeling  of  Cells  Does  not  Significantly  Alter  Cell  Motility 

We  compared  the  migration  speeds  of  wild  type  MDA-MB-231  cells  to  MDA-MB-231 
cells  expressing  a  nuclear  fluorescence  marker  in  the  presence  and  absence  of  EGF  in 
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205  order  to  determine  if  the  introduction  of  nuclear  marker  significantly  impacted  ligand  induced 

206  migration.  Unlabeled  wild  type  cells  were  manual  segmented  and  analyzed  using  Pathfinder, 

207  while  fluorescent  images  were  automatically  segmented  and  analyzed  using  Pathfinder. 

208  Introduction  of  a  fluorescent  nuclear  marker  into  these  cells  did  not  significantly  alter  the 

209  EGF  induced  cellular  migration  speed  (Fig.  S6).  However,  we  do  not  rule  out  that  different 

210  methods  of  gene  delivery  and  types  of  nuclear  markers  (for  instance  a  fluorescent  protein 

211  other  than  histone  H2B)  could  lead  to  permanent  changes  in  cell  migration. 

212 

213 

214  Using  the  average  absolute  angle  of  deflection  to  measure  cellular  persistence 

215  In  order  to  provide  a  means  for  high  throughput  calculation  of  cellular  migration 

216  persistence,  we  used  a  non-traditional,  but  direct,  approach  of  calculating  the  angle  of 

217  deflection  for  each  cell  at  each  time.  Fig.  2B,  bottom  left,  illustrates  how  the  angle  of 

218  deflection  measures  migration  persistence.  The  diagram  represents  a  single  cell,  whose 

219  position  is  measured  at  three  successive  time  points  (1 , 2,  and  3,  respectively).  As  the  cell 

220  travels  from  1  to  2  it  maps  out  a  line  representing  the  trajectory  of  the  cell  between  these 

221  two  times.  Similarly,  as  the  cell  travels  from  2  to  3,  another  line  is  formed.  The  angle  of 

222  deflection  is  the  angle  between  these  two  lines,  where  a  clockwise  turn  has  a  positive  value, 

223  and  a  counterclockwise  turn  has  a  negative  value.  Using  this  calculation,  each  cell  at  each 

224  time  can  be  assigned  an  angle  of  deflection,  such  that  the  sampling  of  many  cells  at  a  single 

225  time  point  can  provide  an  accurate  measurement  of  how  straight  cells  are  migrating  within 

226  the  population.  A  decrease  in  the  average  absolute  value  of  the  angle  of  deflection  for  a 

227  population  of  cells  {(\dDefiection\))  reflects  an  increase  in  the  migration  persistence.  We  use 

228  the  absolute  value  of  the  angle  of  deflection  for  describing  the  persistence  of  migration  in 
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229  large  populations,  rather  than  maintaining  the  sign  of  the  angle  of  deflection,  because  cells 

230  do  not  display  a  bias  in  which  direction  they  prefer  to  turn  (Fig.  3A). 

231 

232  Comparing  methods  to  measure  migration  persistence 

233  Although  we  measure  migration  persistence  using  the  average  absolute  angle  of 

234  deflection,  the  measurement  of  migration  persistence  is  currently  conducted  in  the  cell 

235  migration  field  through  the  determination  of  persistence  time,  which  is  calculated  using  data 

236  fitting  of  time  dependent  mean  squared  displacement  trends  to  Equation  1 ,  where  MSD  is 

237  the  mean  squared  displacement  of  the  cell,  nd  is  the  number  of  dimensions  in  which  cells 

238  are  migrating,  S2  represents  the  squared  speed  of  the  cell,  P  represents  the  persistence 

239  time  of  a  cell,  and  t  is  the  time. 

240 

241  MSD  =  ndS2Pt  [l  -  £  (1  -  e“tp)j  Equation  1 

242 

243  Persistence  time  is  used  to  measure  migration  persistence  because  local  changes  in 

244  mean  squared  displacements  trends  are  likely  to  be  associated  with  changes  in  the  direction 

245  of  cellular  migration,  provided  that  speed  is  taken  into  consideration.  A  key  difference 

246  between  our  approach  and  the  persistence  time  approach  is  that  persistence  time 

247  measurements  focus  on  how  long  a  cell  maintains  a  direction,  while  the  average  absolute 

248  angle  of  deflection  measurements  focus  on  the  degree  to  which  cells  in  a  population  turn  in 

249  each  frame.  Thus,  persistence  time  measurements  reflect  behavior  over  a  specified  interval 

250  of  time  (usually  2-4  hours),  while  the  average  absolute  angle  of  deflection  measurements 

251  reflect  the  relatively  instantaneous  behavior  of  cells. 

252  We  compared  our  technique  of  measuring  migration  persistence  to  the  method  of 

253  measuring  persistence  time  by  examining  the  time-lapse  microscopy  videos  of  two  cell  lines, 
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254  MDA-MB-231  and  HaCaT,  stably  expressing  fluorescent  nuclear  markers  and  treated  with 

255  no  ligand,  Transforming  Growth  Factor  Beta  (TGFP)  or  Epidermal  Growth  Factor  (EGF). 

256  When  such  videos  of  HaCaT  cells  (24-26  hours  post  ligand  stimulation)  are  analyzed  by  the 

257  two  methods,  both  techniques  lead  to  the  same  conclusions;  TGFp  and  EGF  stimulation 

258  cause  increased  migration  persistence,  where  EGF  has  an  impact  of  higher  magnitude  than 

259  that  of  TGFp  (Fig.  3B).  Both  methods  also  agree  when  the  same  analysis  is  applied  to 

260  MDA-MB-231  cells,  where  only  TGFp  has  a  low  magnitude  effect  on  migration  persistence 

261  (Fig.  3B).  The  measured  persistence  time  of  approximately  30  minutes  observed  for  EGF 

262  stimulated  MDA-MB-231  cells  is  consistent  with  similar  results  from  other  studies  [82].  Since 

263  these  two  techniques  yield  the  same  results  under  these  experimental  conditions,  we 

264  conclude  that  both  techniques  accurately  measure  migration  persistence  in  motile  cells. 

265  However,  there  is  one  critical  difference  between  our  method  for  measuring  migration 

266  persistence  compared  to  the  persistence  time  method.  Measuring  the  average  absolute 

267  angle  of  deflection  can  yield  a  measurement  for  migration  direction  for  each  frame,  while 

268  persistence  time  calculations  require  enough  frames  to  construct  a  MSD  vs  time  plot  to  fit  to 

269  Equation  1 .  Thus,  the  persistence  time  calculation,  which  is  traditionally  calculated  from 

270  sampling  over  2-4  hours,  cannot  accurately  report  when  a  cell  turns  and  how  much  it  turns, 

271  but  instead  reports  its  average  tendency  to  turn.  The  value  of  having  the  ability  to  measure 

272  exactly  when  an  to  what  degree  a  cell  turns  will  prove  useful  for  future  investigations  into  the 

273  mechanism  by  which  cells  are  guided  and  will  aid  researchers  in  answering  how  exactly  a 

274  cell  determines  where  to  extend  its  front  end. 

275 

276  Overlapping  intervals  suppresses  noise  in  directional  cellular  behavior 

277  The  Pathfinder  software  is  unique  from  other  software  in  that  it  calculates  angular 

278  information  about  individual  cells.  However,  such  angular  information  requires  that  the 

279  determination  of  a  cell’s  position  be  relatively  noise  free.  When  we  closely  examine  the 
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280  tracks  of  individual  cells,  we  find  that  tracks  exhibit  slight  vibration  on  the  short  timescale  (7 

281  minutes),  such  that  a  cell  that  migrates  relatively  straight  does  not  display  a  perfectly  straight 

282  track.  As  a  result,  we  use  overlapping  intervals  for  our  calculations  of  angular  information  in 

283  order  to  suppress  the  effects  of  such  vibration  on  angular  calculations.  The  schematic 

284  diagram  in  Fig.  4A  illustrates  how  overlapping  intervals  aid  in  the  reduction  of  noise  in  the 

285  calculation  of  cellular  speed,  direction,  and  persistence.  Presented  is  the  path  of  a  single 

286  hypothetical  cell  that  travels  from  positions  1  to  6.  When  calculating  the  trajectory  of  the 

287  movement  from  position  1  to  2,  the  resulting  vector  does  not  accurately  represent  the 

288  underlying  trajectory  of  the  cell  over  time.  However,  as  the  calculation  is  repeated  in  the 

289  same  manner  for  a  change  in  position  from  1  to  increasing  successive  positions,  the 

290  resulting  vectors  quickly  to  converge  on  the  underlying  trajectory  of  the  cell.  Each  interval 

291  represents  the  cellular  behavior  in  a  video  that  has  a  frame  rate  that  is  the  (acquisition  frame 

292  rate)  x  (the  interval  size).  For  instance,  if  the  interval  size  is  3  frames,  then  calculations  are 

293  conducted  on  frames  1 ,4,7,1 0,  and  so  forth.  When  successive  intervals  are  combined,  with 

294  a  single  frame  shift  from  one  interval  to  the  next,  the  resulting  data  provides  a  time 

295  dependent  parameter  that  has  greatly  suppressed  noise. 

296  Fig.  4B  and  C  show  the  effect  of  increasing  interval  size  on  the  time  dependent  average 

297  absolute  angle  of  deflection  trend  and  the  average  standard  deviation  of  the  absolute  angle 

298  of  deflection  trend  for  HaCaT  cells  treated  with  TGFp.  With  an  interval  size  of  1  frame,  both 

299  the  average  (Fig.  4B)  and  the  error  (Fig.  4C)  of  the  angle  deflection  measurements  are 

300  extremely  noisy.  With  increasing  interval  size,  such  noise  is  suppressed,  where  an  interval 

301  size  of  greater  than  2  (corresponding  to  14  minutes)  does  not  yield  significant  additional 

302  suppression  of  noise.  For  all  cellular  experiments  detailed  in  this  investigation,  an  interval 

303  size  of  3  frames  was  used.  This  method  of  overlapping  intervals  was  applied  to  all 

304  measurements,  with  the  exception  of  persistence  time  measurements. 

305 
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306 


307 

308  Measuring  time  dependent  changes  in  migration  persistence  and  speed 

309  Upon  mere  qualitative  assessment  of  cellular  tracks,  cellular  behavior  is  difficult  to 

310  deduce  for  large  populations  of  cells.  For  example,  when  we  examine  MDA-MB-231  and 

311  HaCaT  cell  migration  in  response  to  either  TGFp  or  EGF  by  looking  at  the  tracks  of  cells 

312  between  0  and  35  hours  post  ligand  stimulation,  MDA-MB-231  cells  appear  to  not  change 

313  behavior  in  response  to  ligand  treatments  and  FlaCaT  cells  appear  to  respond  in  the  same 

314  manner  upon  TGFp  and  EGF  stimulation  (Fig.  5A  and  B).  Flowever,  MDA-MB-231  cells  do 

315  in  fact  respond  to  ligand  stimulations  and  HaCaT  cells  do  in  fact  respond  differently  to  TGFp 

316  and  EGF  stimulation  (Video  S7  and  S8),  which  is  elaborated  upon  below.  Thus,  only 

317  through  rigorous  quantitation  can  large  populations  be  characterized  for  their  cellular 

318  behavior. 

319  Using  our  quantitative  approach  to  measuring  cellular  migration,  we  were  able  to 

320  determine  that  either  TGFp  or  EGF  causes  MDA-MB-231  cells  to  migrate  faster,  but  has 

321  almost  no  effect  on  how  persistently  cells  migrate.  Only  stimulation  with  TGFp  causes  a 

322  statistically  significant  (p=  0.003),  but  extremely  small  in  magnitude,  decrease  in  migration 

323  persistence,  as  evident  by  a  slight  elevation  in  the  average  absolute  angle  of  deflection  of 

324  cells  (Fig.  5C,  top  left).  Both  TGFp  and  EGF  treatments  yield  an  increase  in  the  average 

325  speed  of  these  cells,  where  EGF  response  is  early  (approximately  1  hour)  and  TGFp 

326  response  is  late  (approximately  10  hours)  (Fig.  5C,  bottom).  In  contrast  to  MDA-MB-231 

327  cells,  both  speed  and  migration  persistence  are  activated  by  TGFp  and  EGF  treatment  in 

328  HaCaT  cells,  where  the  effects  of  EGF  appear  early  (approximately  1  hour)  and  the  effects 

329  of  TGFp  appear  late  (approximately  1 0  hours)  (Fig.  5D).  Through  our  rigorous  quantitation 

330  method,  we  conclude  that  motility  promoting  ligand  stimulations  can  differ  greatly  from  each 

331  other  in  terms  of  both  the  effect  on  cellular  migration  parameters  and  the  kinetics  of 
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332  activation,  both  of  which  can  be  determined  using  the  Pathfinder  program.  Taken  together, 

333  this  data  illustrates  the  importance  of  time  resolution  in  measuring  cell  migration  responses 

334  to  environmental  stimuli,  as  well  as  the  importance  of  measuring  both  the  speed  and  the 

335  persistence  of  cells  in  order  to  characterize  the  full  behavior  of  cell  migration. 

336 

337  Quantifying  collective  migration  with  nearest  neighbor  analyses  of  angular 

338  measurements 

339  The  increasing  interest  in  the  collective  behavior  of  cells  in  the  cell  migration  field 

340  prompted  us  to  integrate  the  capability  of  characterizing  collective  migration  into  the 

341  Pathfinder  software.  One  already  implemented  method  for  measuring  collective  migration 

342  involves  the  calculation  of  the  size  of  collectively  migrating  streams  of  cells  [33].  Although 

343  this  method  is  certainly  a  valid  way  to  measure  the  ‘collectiveness’  of  cells  within  a 

344  population  of  cells,  it  is  requires  a  predetermined  threshold  to  identify  whether  or  not 

345  neighboring  cells  are  part  of  a  migrating  stream  (a  10  degree  difference  in  migration 

346  direction  or  less  is  required  to  group  cells  together  as  a  collectively  migrating  stream).  In 

347  contrast,  we  used  Pathfinder  and  a  matrix  based  calculation  to  provide  an  alternative 

348  method  for  measuring  how  similar  pairs  of  neighboring  cells  are  in  the  directions  that  they 

349  migrate.  This  method  does  not  require  any  predetermined  thresholds  as  mention  above. 

350  For  each  cell,  Pathfinder  reports  the  angle  of  trajectory,  or  migration  direction.  For  this 

351  parameter,  each  cell  is  assigned  an  angle  that  ranges  from  0  to  359  degrees  relative  to  a 

352  well-defined  set  of  axes  in  the  field  of  view.  Fig.  2B,  right,  illustrates  the  orientation  of  these 

353  axes  in  the  field  of  view,  and  provides  an  example  of  a  cell  that  is  migrating  in  the  direction 

354  of  45  degrees.  In  a  similar  manner  to  the  case  of  angle  of  deflection,  each  cell  gets  such  an 

355  assignment  for  each  time,  providing  time  resolution  of  cellular  direction.  In  order  to 

356  demonstrate  how  the  angle  of  trajectory  calculation  can  be  used  to  characterize  collective 

357  migration,  we  examined  PlaCaT  and  MDA-MB-231  cells  in  confluent  monolayers  in  response 
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358  to  EGF  stimulation.  Manual  inspection  of  the  cellular  tracks  of  these  monolayers  suggests 

359  that  neighboring  cells  migrate  in  a  similar  direction  in  a  ligand  dependent  manner  in  HaCaT 

360  cells  (Fig.  6A,  top).  This  behavior  is  not  qualitatively  observed  for  MDA-MB-231  cells  (Fig. 

361  6A,  bottom).  Since  collective  migration  is  defined  as  the  ability  of  cells  to  adopt  a  common 

362  migration  direction,  we  quantified  the  average  standard  deviation  of  the  angle  of  trajectory, 

363  also  referred  to  as  the  “paired  random  migration  index”  (PRMI  ©Trajectory),  amongst  pairs  of 

364  nearest  neighboring  cells  at  22-26  hours  post  ligand  stimulation.  An  increase  in  this  quantity 

365  indicates  that  nearest  neighboring  cells  are  migrating  in  increasingly  different  directions, 

366  meaning  collective  migration  is  decreasing.  We  excluded  pairs  of  neighbors  in  which  one 

367  cell  migrates  with  a  direction  of  0-90  degrees  and  another  cell  migrates  in  a  direction  of  270- 

368  360  degrees,  as  these  pairs  would  have  a  falsely  high  standard  deviation  due  to  the 

369  discontinuous  transition  between  360  degrees  and  1  degrees.  In  agreement  with  our 

370  qualitative  observations  of  cellular  tracks,  nearest  neighboring  HaCaT  cells  display  a  lower 

371  PRMI  ©Trajectory  in  response  to  EGF  stimulation  (Fig.  6B).  We  compared  nearest  neighbor 

372  behavior  to  random  pairing  behavior  in  order  to  identify  when  shared  behavior  of  cells  is 

373  global  versus  local.  If  random  pairing  does  not  change  the  magnitude  of  the  PRMI  ©Trajectory, 

374  then  the  shared  behavior  is  entirely  global.  In  contrast,  when  such  magnitudes  are  affected 

375  by  random  pairing,  the  phenomenon  is  local.  Upon  random  pairing  of  HaCaT  cells,  a  similar 

376  trend  is  observed  for  the  PRMI  ©Trajectory,  but  the  magnitudes  increase  for  both  mock 

377  treatment  and  EGF  treatment.  Thus,  EGF  stimulation  of  HaCaT  cells  activates  local 

378  collective  migration,  which  diminishes  with  increasing  distance  between  neighbors.  As  a 

379  control,  we  repeated  our  nearest  neighbor  calculations  after  substitution  of  random  angles  of 

380  trajectory  (nearest  neighbors  random  angles)  in  order  to  determine  that  the  maximum  value 

381  for  the  PRMI  ©Trajectory  is  approximately  85  degrees  in  both  the  presence  and  absence  of 

382  EGF  stimulation.  This  maximum  value  is  greater  than  that  of  the  PRMI  ©Trajectory  in  the 

383  absence  of  ligand  stimulation,  revealing  that  HaCaT  cells  do  display  a  small  degree  of 
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384  collective  migration,  which  we  were  not  able  to  detect  upon  manual  qualitative  inspection  of 

385  time-lapse  videos.  Using  the  same  technique  on  MDA-MB-231  cells,  we  found  that  these 

386  cells  display  a  statistically  significant,  but  low  magnitude,  increase  in  the  PRMI  ©Trajectory 

387  amongst  nearest  neighbors  in  response  to  EGF,  suggesting  that  ligand  stimulation  of  these 

388  cells  causes  neighboring  cells  to  exhibit  slight  repulsion,  and  migrate  more  in  opposing 

389  directions  upon  EGF  stimulation  (Fig.  6C).  Random  pairing  of  MDA-MB-231  cells  led  to  an 

390  increase  in  the  magnitude  of  the  PRMI  ©Trajectory,  revealing  that  the  collective  migration  of 

391  these  cells  is  entirely  a  local  phenomenon.  Substitution  of  random  angles  of  trajectories  into 

392  data  sets  revealed  a  similar  maximum  value  for  the  PRMI  ©Trajectory ,  which  was 

393  approximately  85  degrees.  Whether  or  not  the  apparent  collective  migration  behavior  of 

394  neighboring  MDA-MB-231  cells  in  the  absence  of  ligand  stimulation  constitutes  collective 

395  migration  according  to  the  accepted  definition  will  require  further  investigation  into  the 

396  requirement  of  cellular  junctions  in  this  process.  Flowever,  it  is  worth  noting  that  the 

397  neighboring  cells  do  in  fact  have  physical  contact  with  each  other  (Fig.  S9),  albeit  for  short 

398  timescales  (data  not  shown). 

399 

400  Using  nearest  neighbor  analyses  of  angular  measurements  to  characterize  collective 

401  migration  in  epithelial  sheets 

402  In  order  to  determine  the  versatility  of  our  method  of  collective  migration  quantification, 

403  we  repeated  our  experiments  on  FlaCaT  cells  using  EGF  stimulation,  but  under  the  condition 

404  in  which  cells  were  arranged  into  epithelial  sheets  (Fig.  6D  and  Video  S10).  We  asked 

405  whether  or  not  ligand  stimulation  instantaneously  affects  all  cells  equally  in  an  epithelial 

406  sheet,  and  found  that  the  collective  behavior  response  to  EGF  stimulation  propagates  away 

407  from  the  leading  edge  of  an  epithelial  sheet  overtime  as  indicated  by  propagation  of 

408  decreasing  PRMI  ©Trajectory  away  from  the  leading  edge  of  an  epithelial  sheet  (Fig.  6E).  Thus, 
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409  our  technique  of  quantitatively  characterizing  collective  migration  has  proven  useful  to 

410  determine  the  spatial  collective  migration  properties  throughout  a  population  of  cells. 

411 

412 

413 

414  Discussion 

415 

416  The  pathfinder  software  provides  novel  useful  tools  for  studying  individual  and 

417  collective  cell  migration 

418  Our  Pathfinder  software  is  the  first  high  throughput  automated  cell  migration  software  to 

419  conduct  measurements  of  instantaneous  angular  parameters  for  the  characterization  of  cell 

420  migration  behavior.  Our  parameter  of  the  angle  of  deflection  has  proven  useful  for 

421  characterizing  the  migration  persistence  of  cells  in  two  dimensional  fluorescence  microscopy 

422  experiments.  In  addition  to  helping  researchers  characterize  the  molecular  mechanisms 

423  behind  migration  persistence,  this  parameter  can  be  used  to  identify  the  differential  temporal 

424  responses  of  cells  to  distinct  ligand  stimulations,  and  can  provide  insight  into  the  molecular 

425  mechanisms  behind  each  ligand  response.  Furthermore,  the  available  cell  migration 

426  software  has  been  entirely  devoid  of  tools  that  can  be  used  to  readily  characterize  collective 

427  migration.  As  a  result,  investigations  of  collective  migration  have  been  more  appropriately 

428  described  as  qualitative  rather  than  quantitative.  The  Pathfinder  software  provides  the  first 

429  high  throughput  tool  to  quantitatively  characterize  collective  cell  migration  using  the  angle  of 

430  trajectory  parameter,  which  can  be  used  to  determine  the  “Paired  Random  Migration  Index” 

431  (PRMI  ©Trajectory)  amongst  nearest  neighbors.  This  tool  is  likely  to  propel  the  collective 

432  migration  field  forward,  as  it  can  provide  a  means  for  accurately  measuring  the  degree  to 

433  which  cells  are  migrating  collectively,  with  both  spatial  and  temporal  resolution.  Our  method 
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434  of  measuring  collective  migration  was  able  to  reveal  hidden  collective  cellular  behavior  that 

435  cannot  be  easily  detected  qualitatively  from  cellular  tracks,  such  as  the  local  collective 

436  migration  behavior  amongst  pairs  of  MDA-MB-231  cells,  and  the  ability  of  EGF  stimulation  to 

437  suppress  this  behavior.  In  conclusion,  our  Pathfinder  software  provides  novel  techniques  for 

438  characterizing  cellular  migration  in  a  high  throughput  platform,  which  are  likely  to  aid  the  cell 

439  migration  field  in  its  investigation  of  the  molecular  mechanisms  behind  individual  and 

440  collective  migration. 

441 
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457  Figure  1.  The  Pathfinder  cell  motility  program  uniquely  incorporates  measurements  of 

458  cellular  position,  speed,  direction,  and  persistence. 
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459  Although  several  cellular  motility  programs  are  available  to  measure  cellular  motility  in  terms  of 

460  position  and  speed,  only  the  Pathfinder  program  additionally  reports  both  cellular  migration 

461  direction  and  cellular  migration  persistence. 

462 

463  Figure  2.  Angular  measurements  of  cellular  migration  can  reveal  cellular  behavior. 

464  A)  The  Pathfinder  program  converts  time-lapse  microscopy  videos  of  fluorescent  HaCaT  cells 

465  (left)  to  cellular  tracks  (right).  B)  The  Pathfinder  program  uses  the  positional  information  of 

466  cellular  tracks  to  calculate  speed,  where  dS  represents  change  in  position  and  eft  represents 

467  change  in  time  (top),  migration  persistence  through  the  absolute  angle  of  deflection  (bottom), 

468  and  the  migration  direction  relative  to  a  well-defined  axis  orientation  in  the  field  of  view  (right). 

469 

470  Figure  3.  Average  absolute  angle  of  deflection  measurements  accurately  depict  the 

471  migration  persistence  of  cells. 

472  A)  Cells  do  not  prefer  to  turn  right  or  left  in  either  the  presence  (right)  or  absence  (left)  of  EGF 

473  stimulation.  A  binned  histogram  of  percent  of  cells  versus  percent  of  right  turns  is  normal  and 

474  centered  around  50  percent.  B)  A  comparison  of  persistence  time  calculations  and  average 

475  absolute  angle  of  deflection  {(\0 Deflection]))  methods  for  measuring  migration  persistence  yields 

476  identical  trends  in  both  the  presence  and  absence  of  either  TGFfB  or  EGF  for  MDA-MB-231  cells 

477  and  PlaCaT  cells.  Double  asterisks  indicate  a  p  value  <  0.01.  Each  condition  represents  greater 

478  than  200  cells  for  persistence  time  measurements,  and  greater  than  1 000  cells  for  average 

479  absolute  angle  of  deflection  measurements. 

480  Figure  4.  Overlapping  Intervals  Suppress  Angular  Noise  in  Cellular  Migration. 

481  A)  A  schematic  representation  of  a  cellular  track  illustrates  how  increasing  interval  size  results 

482  in  the  calculation  of  a  cellular  migration  direction  that  is  resistant  to  track  vibration  noise.  B)  The 

483  average  absolute  angle  of  deflection  as  a  function  of  time  for  TGFp  treated  FlaCaT  F12B 

484  mCherry  cells  for  an  interval  size  of  one  frame  shows  strong  scattering  of  measurements.  Upon 

485  increasing  interval  size,  such  scattering  is  suppressed.  C)  Such  scattering  can  be  also 

486  measured  by  measuring  the  standard  deviation  in  the  absolute  angle  of  deflection,  which  can  be 

487  suppressed  in  a  similar  manner  by  increasing  interval  size.  Data  represents  greater  than  1000 

488  cells  for  each  plot. 

489 

490 

491 

492  Figure  5.  Measuring  individual  cellular  behavior  with  speed  and  migration  persistence 

493  reveals  cell  type  and  ligand  specific  cellular  migration  behavior. 
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494  Cellular  tracks  of  low  density  cells  are  displayed  for  treatments  of  either  Mock,  TGFp  or  EGF  for 

495  MDA-MB-231  cells  (A)  or  HaCaT  cells  (B).  Calibration  bars  represent  150  pm.  C)  Neither  TGF- 

496  Beta  nor  EGF  stimulation  affects  migration  persistence  in  MDA-MB-231  cells  (top).  In  contrast, 

497  both  treatments  affect  cellular  speed,  but  with  different  induction  kinetics  (bottom).  D)  In  HaCaT 

498  cells,  both  ligand  treatments  affect  migration  persistence  and  cellular  speed  (top  and  bottom, 

499  respectively).  However,  EGF  stimulates  migration  persistence  with  earlier  kinetics  than  that  of 

500  TGFp  (top),  and  EGF  is  a  poor  stimulator  of  migration  speed  (bottom,  right).  Each  condition 

501  represents  greater  than  1000  cells. 

502 

503  Figure  6.  Angular  measurements  can  be  used  to  quantify  collective  migration  behavior. 

504  A)  Cellular  tracks  of  confluent  monolayers  of  HaCaT  (top)  and  MDA-MB-231  cells  (bottom)  in 

505  the  presence  and  absence  of  EGF  stimulation.  Calibration  bar  represents  150  pm.  B)  Confluent 

506  monolayers  of  HaCaT  cells  in  the  presence  and  absence  of  EGF  stimulation  were  quantified  for 

507  their  collective  migration  behavior  by  calculating  the  average  standard  deviation  of  the  angle  of 

508  trajectory  (also  called  the  paired  random  migration  index  (PRMI  0Tajectory)  amongst  nearest 

509  neighboring  cells.  Random  pairing  was  used  to  determine  whether  the  observed  behavior  was 

510  local  or  global  amongst  the  population.  C)  The  same  quantification  was  conducted  for  MDA-MB- 

511  231  cells.  D)  Cellular  tracks  of  epithelial  sheets  of  HaCaT  cells  in  the  presence  and  absence  of 

512  EGF  stimulation.  Calibration  bar  represents  150  pm.  E)  Inspection  of  the  spatial  distribution  of 

513  collective  migration  behavior  reveals  that  EGF  stimulation  elicits  collective  migration  that 

514  propagates  away  from  the  leading  edge.  Double  asterisks  indicate  a  p  value  <  0.01 . 

515 

516  Supplementary  Data  SI.  Pathfinder  Program. 

517  The  Pathfinder  program  is  an  executable  JAVA  program  that  requires  installation  of  JAVA 

518  Runtime  Environment  on  a  Windows  operating  system. 

519 

520 

521  Supplementary  Data  S2.  Parameter  descriptions  for  the  Pathfinder  program  GUI. 

522  User  input  parameters  are:  Cell  Outer  Radius,  Cell  Minimum  Radius  (Cutoff),  Percentage  of 

523  Pixels  with  Nuclear  Signal  (Percentile),  How  Far  an  Average  Cell  is  Tolerated  to  Migrate  From 

524  Frame  to  Frame  (Disp./Frame),  How  Many  Frames  to  Bin  for  Calculations  (Frame  Binning), 

525  Minimum  Tack  Length  (Min.  Trajectory  Length),  Number  of  Parallel  Threads  (Threads),  Folder 

526  Path  for  Folder  with  Videos  (AVI  Folder). 

527  Supplementary  Data  S3.  A  description  of  output  calculations  from  the  Pathfinder 

528  program. 
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529  Each  cell  receives  a  cellular  ID  number  (1),  for  each  frame  (2).  In  each  frame  a  cell  is  assigned 

530  an  X  (3)  and  Y  (4)  coordinate,  a  displacement  from  the  last  frame  in  pixels  (5),  an  angle  of 

531  trajectory  (6),  an  angle  of  deflection  (7)  and  a  mean  squared  displacement  (8).  Mean  squared 

532  displacements  can  be  used  to  calculate  the  persistence  time  for  a  cell.  For  the  population  of 

533  cells,  Pathfinder  reports  the  frame  (9)  dependent  change  in  the  average  displacement  (10),  the 

534  average  angle  of  trajectory  (11),  the  percentage  of  cells  turning  greater  than  90  degrees  (12), 

535  and  the  average  absolute  angle  of  deflection  (13).  Additionally,  Pathfinder  reports  a  binned 

536  histogram  of  percent  of  cells  versus  the  possible  migration  directions  from  0  to  359  degrees  (14 

537  and  1 5).  Lastly,  the  number  of  cellular  tracks  is  reported  (16). 

538  Supplementary  Data  S4.  A  sample  output  from  Pathfinder. 

539  A  sheet  of  PlaCaT  H2B-mCherry  cells  were  analyzed  for  individual  cell  migration  from  a 

540  microscopy  video  between  22  and  25.5  hours  post  100  nM  EGF  stimulation.  Frames  were 

541  acquired  every  7  minutes.  Analysis  was  conducted  with  a  frame  binning  of  3  frames  (21 

542  minutes). 

543  Supplementary  Data  S5.  A  MATLAB  script  to  convert  mean  squared  displacement  (MSD) 

544  versus  time  data  from  a  Pathfinder  output  into  persistence  time. 

545  This  MATLAB  script  requires  user  input  of  MSD  and  time  data  for  each  cell  and  converts  this 

546  information  into  persistence  time. 

547  Supplementary  Data  S6.  Wild  type  MDA-MB-231  cells  and  MDA-MB-231  H2B-mCherry 

548  cells  migrate  with  similar  speeds  in  the  presence  and  absence  of  EGF  stimulation. 

549  Brightfield  microscopy  videos  of  mock  and  EGF  treated  wild  type  (WT)  MDA-MB-231  cells  were 

550  manually  measured  for  position  over  the  course  of  a  1 0  frame  interval  (7  minutes/frame)  after  24 

551  hours  ligand  or  mock  stimulation  and  the  average  speed  of  cells  was  calculated  with  a  frame 

552  binning  of  3.  The  same  analysis  was  done  on  parallel  videos  of  MDA-MB-231  using  pathfinder, 

553  which  yielded  similar  results  for  WT  and  labeled  cells  in  the  speed  of  migration  in  the  presence 

554  and  absence  of  EGF.  50  cells  were  used  for  this  comparison  for  each  condition. 

555  Supplementary  Video  S7.  MDA-MB-231  cells  at  low  density  upon  either  mock,  TGFp,  or 

556  EGF  treatment. 

557  MDA-MB-231  cells  with  an  FI2B-mCherry  nuclear  marker  were  observed  by  time-lapse 

558  microscopy  using  the  mCherry  fluorescence  channel.  Each  frame  represents  7  minutes. 

559  Supplementary  Video  S8.  HaCaT  cells  at  low  density  upon  either  mock,  TGFp,  or  EGF 

560  treatment. 

561  FlaCaT  cells  with  an  FI2B-mCherry  nuclear  marker  were  observed  by  time-lapse  microscopy 

562  using  the  mCherry  fluorescence  channel.  Each  frame  represents  7  minutes. 

563 
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564  Supplementary  Data  S9.  MDA-MB-231  cells  maintain  physical  contact  with  their  nearest 

565  neighboring  cell. 

566  Brightfield  microscopy  of  EGF  treated  MDA-MB-231  cells  reveals  that  nearest  neighboring  cells 

567  have  physical  contact  with  each  other. 

568 

569  Supplementary  Video  S10.  Epithelial  sheets  of  HaCaT  cells  upon  either  mock  or  EGF 

570  treatment. 
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HaCaT  cells  with  an  H2B-mCherry  nuclear  marker  were  assembled  into  epithelial  sheets  and 
observed  by  time-lapse  microscopy  using  the  mCherry  fluorescence  channel.  Each  frame 
represents  7  minutes. 
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